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and  coded  the  models  used  to  predict  the  nuclear  thermal  radiation 
environments.  Mr.  John  Dishon  provided  the  description  and  data  on  the 
SAI  thermal  radiation  simulator,  as  well  as  fruitful  discussions  on 
results  of  past  experiments  he  performed. 

Special  thanks  are  given  to  Dr.  Walter  F.  Dudziak  of  Information 
Science,  Inc.,  for  providing  thermal  radiation  data  from  some  nuclear 
shots.  The  authors  also  thank  Mr.  Tom  Kennedy,  who  monitored  the  tech¬ 
nical  performance  of  this  effort  for  the  Defense  Nuclear  Agency,  for 
his  helpful  guidance  and  suggestions. 
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Conversion  factors  for  U.S.  customary 
to  metric  (SI)  units  of  measurement. 


To  Convert  From 

To 

Multiply  By 

atmosphere  (normal 

kilo  pascal  (kPa) 

'  1.013  25  X  E  +2 

bar 

kilo  pascal  (kPa) 

1.000  000  X  E  +2 

calorie  (thermochemical) 

joule  (J) 

4.184  000 

o 

cal  (thermochemical  )/cm 

?  ? 
mega  joule/m  (MJ/m  ) 

4.184  000  X  E  -2 

erg 

joule  (J) 

1.000  000  X  E  -7 

erg/second 

watt  (W) 

1 .000  000  X  E  -7 

foot 

meter  (m) 

3.048  000  X  E  -1 

inch 

meter  (m) 

2.540  000  X  E  -2 

kilotons 

terajoul es 

4.183 

pound-force/inch^  (psi) 

kilo  pascal  (kPa) 

6.894  757 
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SECTION  1 


INTRODUCTION 

An  important  difference  between  a  nuclear  explosion  (NE)  and  a 
chemical  high  explosive  (HE)  is  the  relatively  large  amount  of  energy 
released  during  the  nuclear  explosion  in  the  form  of  thermal  radiation. 
Synergistic  effects  due  to  air  blast  and  thermal  radiation  are  known  to 
exist.  These  effects  have  been  difficult  to  obtain  in  the  field  pri¬ 
marily  due  to  the  lack  of  large  scale  thermal  simulation  devices  that 
produce  significant  radiative  output. 

SAI  has  been  developing  a  large  scale  thermal  radiation  simu¬ 
lator  for  DNA,  and  it  seems  feasible  to  now  consider  fielding  such 
arrays  on  a  planned  HE  event,  MISTY  CASTLE.  This  program  was  conducted 
to  develop  data  and  guidelines  that  will  assist  future  experimental 
planning. 

SAI  accomplished  two  tasks  in  preparation  for  planning  an  HE 
field  simulation  of  a  low  yield  nuclear  event  detonating  on  the  ground. 
The  objective  of  the  first  task  was  to  provide  a  description  of  the 
thermal  radiation  environments  that  are  to  be  simulated  in  the  HE  event. 
The  objective  of  the  second  task  was  to  provide  a  first  look  at  how  an 
ensemble  of  DNA  thermal  simulation  devices  might  be  positioned  to 
achieve  an  acceptable  simulation  at  various  ground  ranges  from  the  HE 
detonation  site  without  seriously  distorting  the  airblast  flow  fields. 
The  results  from  these  two  tasks  are  expected  to  help  guide  subsequent 
test  planning  and  lay  the  foundation  for  estimating  reasonable  expecta¬ 
tions  for  success  in  achieving  combined  thermal  and  airblast  environ¬ 
ments  . 
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SECTION  2 


THERMAL  RADIATION  ENVIRONMENT  FROM  A  NUCLEAR  1  KT  SURFACE  BURST 

A  general  description  of  the  radiative  output  characteristics 
from  a  nuclear  weapon  can  be  found  in  references  1  and  2.  However, 
more  detailed  theoretical  predictions  of  thermal  radiation  from  free- 
air  nuclear  explosions  have  been  modeled  and  were  used  for  this  effort. 
The  models  for  the  fireball  thermal  radiative  output  are  described  in 
Section  2.1.  How  these  were  utilized  is  described  in  Section  2.2. 

Results  of  calculations  using  these  models  are  presented  in  Section  2.3 
(and  the  appendices).  Finally,  in  Section  2.4  a  discussion  is  presented 
of  possible  uncertainties  that  may  exist  from  effects  not  modeled.  An 
example  of  such  an  effect  is  the  interaction  of  the  radiation  with  the 
ground  and  how  this  might  influence  the  radiative  output  of  the  fireball. 
This  effect  is  important  to  this  effort  since  the  simulation  is  of  a 
surface  burst. 

2.1  THE  AFWL  PREDICTIONS  AND  MODELS  OF  THE  THERMAL  RADIATION  FROM 

A  NUCLEAR  BURST 

The  Air  Force  Weapons  Laboratory  conducted  a  seven  year  theo¬ 
retical  effort  (1963-1970)  to  predict  early  time  fireball  growth* 

(i.e.,  times  from  x-ray  deposition  into  the  cold  air  surrounding  the 
nuclear  device  to  a  time  at  least  ten  times  that  of  second  maximum) 
using  a  code  called  SPUTTER  (see  references  3  and  4).  As  a  result  of 
the  AFWL  effort,  detailed  radiative  and  hydrodynamic  environments  could 
be  predicted  for  singular  nuclear  weapons  bursting  at  various  altitudes 
from  0  to  60  km  in  the  atmosphere.  Weapon  yields  ranged  from  20  tons 
to  100  megatons.  The  predictions  were  free-field  (i.e.,  no  account  was 


*Under  sponsorship  by  DASA,  ARPA  and  the  Air  Force. 


13 


taken  of  presence  of  the  ground)  computer  calculations,  which  now  serve 
as  an  extensive  fireball  data  base.  They  agree  well  with  experimental 
data  taken  in  the  late  50' s  and  early  60 's  in  the  U.S.  above-ground 
nuclear  testing  (references  5  and  6). 

Some  of  the  important  features  of  the  SPUTTER  results  were 
then  modeled  at  the  AFWL  in  the  early  70' s  in  order  to  obtain  simple 
algorithms  (references  7  and  8).  For  example  at  sea  level  the  power 
out  at  second  maximum  from  a  free-field  nuclear  fireball  was  found  to 
be  well  represented  by 

P  =  3.53  x  1012W°*59  cal/sec  (1) 

where  W  =  total  yield  in  kt.  Similarily  the  time  to  second  maximum  at 
sea  level  is  well  represented  by 

^2max  =  seconds.  (2) 

For  1  kt  the  AFWL  algorithm  agrees  quite  well  with  the  SPUTTER  calcula¬ 
tion  at  0  km.  Table  1  lists  the  SPUTTER  calculations  serving  as  a  data 
base  for  the  AFWL  modeling  study  of  radiative  output.  Each  entry  repre¬ 
sents  a  detailed  one-dimensional  radiative  hydrodynami cs  computation 
taking  over  50  hours  on  a  CDC  6600  using  SPUTTER.  The  calculations 
employ  many  electromagnetic  frequency  groups  (typically  20)  to  repre¬ 
sent  thermal  radiation  since  grey  body  approximations  were  found  to  be 
inadequate  (reference  5  ) ;  hence  a  large  computation  time  results  when 
coupled  with  zoning  requirements  for  modeling  the  hydrodynamics. 

The  agreement  of  the  AFWL  algorithm  for  power  and  time  at  second 
maximum  with  SPUTTER  is  quantitatively  good  as  indicated  by  Figures  1 
and  2.  Where  the  SPUTTER  data  are  considered  reliable,  the  AFWL  models 
are  good  representations.  Results  of  the  FB21  SPUTTER  calculation  used 
in  the  AFWL  study,  are  included  and  show  results  for  a  1  kt  free-field 
(no  yield  doubling,  see  Section  2.3)  burst  detonating  at  sea  level 
(0  km  al ti tude) .* 


*These  results  have  to  be  scaled  to  be  used  for  a  1  kt  surface  burst. 
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Table  1.  SPUTTER  runs  used  in  generating 
AFWL  thermal  output  model  . 

Extracted  from  Reference  7 


Ident 
(FB  No.) 

Simulated  Yield 
(kt) 

Detonation  Alt 

Air  Ratio  of  Air  Altitudp 

Density  Density  to  Density  /km\ 
(gm/cm)  at  Sea  Level 

20 

0.1  (0.094)* 

1 . 225E-03 

1 

o** 

21 

1  0) 

1 . 225E-03 

1 

0 

2 

10  (10) 

8. 253E-05 

0.06737 

20.3 

14 

30  (30) 

1 . 225E-03 

1 

0 

11 

30  (30) 

1 .161E-04 

0.09477 

18.3 

23 

30  (30) 

1 . 355E-05 

0.01106 

32 

19 

30  (30) 

1 .351E-06 

0.001103 

48.3 

7A 

198  (198) 

1 . 225E-03 

1 

0 

13 

198  (199) 

4.587E-04 

0.37445 

9.1 

3F 

198  (200) 

1 .308E-04 

0.106776 

18.3 

5 

198  (200) 

1 . 356E-05 

0.01107 

32 

S2 

198  (199) 

7.2579E-06 

0.0059248 

36 

S5 

198  (200) 

4.6267E-06 

0.0037769 

39 

15 

198  (198) 

3 . 996E-06 

0.003262 

40 

S4 

198  (198) 

3.4564E-06 

0.002821 

41 

S3 

198  (197) 

2 . 995E-06 

0.002445 

42 

SI 

198  (197) 

2 . 2589E-06 

0.001844 

43 

IB 

198  (198) 

1 .351E-06 

0.001103 

48 

12 

1000  (993) 

1 .161E-04 

0.09477 

18.3 

10 

3800  (3843) 

1 .225E-03 

1 

0 

18 

3800  (3776) 

1.161E-04 

0.09477 

18.3 

9 

3800  (3773) 

2.780E-06 

0.0022694 

43 

26 

3800  (3773) 

1 .269E-06 

0.001036 

48.3 

25 

155000  (155000) 

1 .161E-04 

0.09477 

18.3 

^Numbers  in  parentheses  are  yields  based  on  the  total  energy  as 
determined  during  the  actual  SPUTTER  CALCULATION. 


**This  altitude  corresponds  to  sea  level.  However,  there  is  no  account 
of  ground  in  the  SPUTTER  calculations,  they  are  free-field. 
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Time  to  Second  Maximum  (Sec) 


Sharp's  Fit  (Sea  Level) 
t  -0. 038W0 • 44  sec 
max  (W  in  Kt) 


Figure  1.  Sharp's  fit  to  SPUTTER  data  compared  with  SPUTTER 
results  (at  sea  level)  for  time  to  second  maximum. 
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Power  Out  (Cal/Sec) 


Sharp's  Fit:  (Sea  Level  Power  Out  3.53E12W0'^ 
cal  /  sec)  (W  in  Kt) 


Figure  2.  Sharp's  fit  to  SPUTTER  data  compared  with  SPUTTER 
results  (at  sea  level)  for  power  out  at  second 
maximum. 
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Figure  3  presents  the  radius  of  various  pressure  fronts  versus 
time  for  this  same  calculation,  and  is  shown  as  a  frame  of  reference  to 
briefly  familiarize  the  reader  with  the  calculation.  The  very  early 
fireball  is  about  4  m  in  radius.  As  the  fireball  develops  and  grows 
radiatively,  the  air  in  front  of  the  hot  radiating  fireball  starts  to 
shock,  and  as  this  air  shock  develops  it  becomes  optically  opaque  so 
that  for  some  time  the  fireball  and  air  shock  appear  to  grow  together. 

At  about  0.07  msec  the  faster  moving  debris  shock  catches  up  with  the 
main  air  shock  and  this  interaction  reinforces  the  main  shock  reflecting 
another  shock  back  towards  the  burst  point.*  At  some  time  after  debris 
shock  arrival  (not  obvious  from  the  figure)  the  air  shock  breaks  away 

from  the  fireball.  Both  continue  to  grow,  the  fireball  slower  than  the 

/ 

air  shock.  Finally,  fireball  growth  is  slowed  appreciably  as  it  comes 
into  pressure  equilibrium  and  the  main  air  shock  continues  on.  At  some 
intervening  time  (tens  of  msec)  the  shock  no  longer  is  optically  opaque 
and  the  true  interior  fireball  can  be  seen.  Until  such  time,  the 
apparent  fireball  was  the  air  shock.  Figure  4  presents  the  correspond¬ 
ing  time  sequence  of  power  out  from  the  FB21  fireball  and  it  clearly 
shows  the  second  maximum.  First  minimum  occurs  near  3  milliseconds. 

Figure  5  is  the  total  energy  radiated  away  from  the  fireball;  approxi¬ 
mately  35%  of  the  yield  has  been  thermally  radiated  by  1  second.  Data 
shown  for  times  earlier  than  1  millisecond  are  unreliable,  as  will  be 
discussed  later. 

The  AFWL  model  also  provides  the  time  and  yield  dependence  of 
the  radiative  output.  Figure  6  shows  the  power  radiated  from  10-^  seconds 
to  100  seconds  for  seven  different  yields.  The  power  radiated  out  is  seen 
to  be  yield  dependent  although  the  ratio  of  total  radiation  to  total  yield 


*This  shock  "reflects"  from  the  center  and  moves  outward  again.  The 
process  repeats,  with  shocks  becoming  weaker,  and  leads  to  many  weak 
pressure  fronts  seen  in  Figure  3.  The  coarse  zoning  tends  to  smear 
these  fronts  leaving  numerical  noise.  The  criteria  for  representing 
a  pressure  front  are  based  on  the  relative  values  of  the  artificial 
viscosity  pressure,  P2,  and  the  equation-of-state  pressure,  PI. 
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Figure  4.  FB21  power  out  (cal/sec). 
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Figure  5.  FB21  energy  radiated. 
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Figure  6.  Thermal  power  from  free-field  burst  at  sea  level 
as  a  function  of  yield.  (Extracted  from 
reference  7). 
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1x10 


is  relatively  the  same  at  both  small  and  large  yields.  However,  even 
though  the  total  radiation  to  yield  ratio  is  fairly  constant,  it  is 
inappropriate  to  conclude  that  total  fluence  will  scale  with  yield  in 
the  same  fashion  as  the  hydrodynamics  scale.  As  can  be  easily  demon¬ 
strated,  higher  yields  will  produce  higher  total  fluence*at  the  same 

1  /3 

hydrodynamical ly  scaled  range;  this  total  fluence  will  scale  as  W 
Consider  a  nuclear  device  of  yield,  W,  detonating  above  the  ground  at 
scaled  height,  h,  and  ground  range,  r.  Distance  in  all  directions  can 
be  scaled,  where 


Scaled  distance 


true  distance 

w1/3 


This  scaling  is  useful  since  airblast  at  any  position  becomes  relatively 

insensitive  to  yield.  But  the  fluence  arriving  at  comparatively  scaled 

positions  is  proportional  to  the  product  of  the  thermal  efficiency,  e, 

2 

and  yield,  W,  divided  by  R  .**  Since  e  is  only  weakly  dependent  on  yield,*** 

1  /3 

total  fluence  remains  approximately  proportional  to  W  at  a  fixed  hydro- 
dynamical  ly  scaled  range. 

Instantaneous  flux  is  not  so  simply  scaled,  since  the  time 
phasing  of  the  radiated  power  is  complex.  If  the  time  to  second  maxi¬ 
mum  is  proportional  to  then  the  time  shift  in  hydrodynamical ly 

scaled  time****will  be  W°  ’  ^  .  The  flux  at  second  maximum  scales  as 
when  comparing  hydrodynamical ly  scaled  ranges.  A  comparison 
can  be  made  between  the  data  in  Figure  4  and  the  1  kiloton  curve  in 
Figure  6  (1  cal/sec  =  4.2  watts)  and  is  shown  as  Figure  7.  The 


*Fluence  is  sometimes  called  radiant  exposure  and  is  the  energy  per 
unit  area  incident  on  some  target  surface.  Conventionally  its  units 
are  cal/cm^. 


**Since  the  height  of  burst  and  ground  range  have  all  been  proportion¬ 
ally  scaled  as  l\P'  ,  the  angle  of  arrival  for  the  radiation  is  scale 
independent. 

***The  AFWL  fit  for  sea  level  efficiency  is  0.3505V/3’*3333  where  W  is 
the  yield  in  kilotons  (reference  7). 


****Scaled  time 


true  time 
W1'3 
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SPUTTER  data  represented  by  the  circles,  are  not  shown  for  first  max¬ 
imum  because  they  are  not  considered  reliable.  Inspection  of  the  very 
early  time  SPUTTER  power  data  by  the  authors  (not  shown  in  this  report) 
shows  that  the  radiated  power  oscillates  rapidly  about  the  expected 
values.  This  is  due  to  the  zoning  used  in  the  SPUTTER  calculations 
which  is  inadequate  to  accurately  represent  radiation  fronts  at  those 
times.  Since  the  mean  free  path  of  the  dominant  radiation  at  early 
times  is  of  the  order  of  tenths  of  centimeters,  zones  of  tens  of  meters 
in  size  will  either  be  too  transparent  or  opaque  and  will  rapidly  change 
from  transparent  to  opaque  to  transparent  again  as  the  radiation  front 
moves  through  any  cell.  This  will  alternately  allow  too  much  radiation 
to  leak  out,  and  then  too  little.  For  the  first  maximum,  nuclear  data 
were  used  to  suggest  the  relationships  established  (see  reference  7). 

The  agreement  between  the  AFWL  model  and  the  detailed  SPUTTER 
power  data  for  FB21  are  reasonable,  although,  there  are  significant 
variations  beyond  the  time  of  second  maximum  as  was  seen  in  Figure  7. 

The  AFWL  model  produces  better  agreement  with  the  SPUTTER  data  for  large 
yields  as  seen  in  Figure  8  showing  a  200  kiloton  calculation.  However, 
this  effort  addresses  low  yields,  so  the  effect  of  these  inaccuracies 
is  considered. 

The  spectral  content  of  the  radiation  is  also  yield  dependent. 
This  is  shown  in  Figure  9,  which  also  shows  the  spectral  content  of 
three  blackbodies  (6000,  3000,  and  2450  degrees  Kelvin).  This  figure 
represents  another  AFWL  model  of  the  SPUTTER  data.  The  SPUTTER  cal¬ 
culations  included  the  calculation  of  the  radiative  transport  through 
the  fireball  for  a  number  of  intervals  in  the  electromagnetic  spectrum 
(typically  20),  hence  spectral  data  are  available  averaged  over  these 
intervals.  Details  of  the  intervals  used  can  be  found  in  references  4, 
5,  6  and  8.  For  comparison  the  same  two  representative  calculations 
shown  above  were  chosen  to  illustrate  the  agreement  between  the  AFWL 
model  of  the  time  dependence  of  the  spectrum.  FB21  (1  kt  at  sea  level) 
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Figure  7.  Thermal  power  for  FB-21  (1  kt  at  sea  level). 
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Figure  8.  Thermal  power  for  FB-7A  (198  kt  at  sea  level) 
(Extracted  from  reference  7). 
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Figure  9.  Comparison  of  function  result  with  blackbody  spectra  for 
0  km  burst  altitude. 

(Extracted  from  reference  8,  Figure  64.) 
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is  shown  as  Figure  10  and  FB7A  (200  kt  at  sea  level)  is  shown  as 
Figure  11.  The  circles  are  again  the  specific  SPUTTER  data  used  by  the 
AFWL  in  the  generation  of  their  models.  Here  again  the  agreement  for 
the  200  kiloton  yield  is  definitely  better  than  that  for  1  kiloton. 

The  SPUTTER  calculation  suggests  a  hotter  fireball  for  1  kiloton  than 
obtained  from  the  AFWL  model. 

Additional  information  on  these  calculations  and  models  can  be 
found  in  References  3  through  8. 

2.2  CALCULATION  OF  THE  THERMAL  RADIATION  ENVIRONMENTS  ON  THE 

GROUND 

Once  the  details  of  the  radiative  output  of  the  nuclear  weapon 
fireball  are  available,  the  problem  of  calculating  the  radiation  reach¬ 
ing  a  detector  located  outside  of  the  fireball  becomes  one  of  appropri¬ 
ately  integrating  over  the  sources  of  radiation,  assuming  progapation 
path  absorption  and  scattering  are  insignificant  (references  9  and  10). 
For  the  ranges  involved  here  (i.e.,  250  m  or  closer)  for  1  kiloton, 
total  power  absorption  is  insignificant  on  a  clear  day  when  visibility 
is  high.  When  visibility  is  low  it  can  become  significant.  For  example, 
when  the  visibility  is  1  mile  the  transmittance  is  about  0.8  for  a  range 
of  250  m. 


For  the  effort  reported  here  the  fireball  was  approximated  as 
an  isotropic  surface  emitter.  Since  the  interest  at  DNA  is  simulating 
a  1  kiloton  surface  burst,  certain  geometric  simplifications  were  pos¬ 
sible  when  integrating  over  the  surface  of  the  fireball  to  obtain 
detector  radiation.  These  simplifications  are  discussed  within  Sec¬ 
tion  2.2.3.  The  following  derivation,  which  is  more  general  than  that 
needed  for  the  surface  burst  case,  is  for  a  detector  whose  normal  is  in 
the  plane  perpendicular  to  the  earth  and  containing  the  points  of  detec¬ 
tion  and  burst. 
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Figure  10.  Thermal  spectrum  for  FB-21  (1  kt  at  sea  level). 
(Extracted  from  reference  8) 
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Fractional  Energy/Micron 


Figure  11.  Thermal  spectrum  for  FB-7A  (198  kt  at  9  km). 
(Extracted  from  reference  8). 
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2.2.1  Formulation  for  the  Truncated  Radiating  Sphere 

A  field  point  P  (Figure  12)  is  positioned  on  a  ground  plane  a 
distance  r  from  the  center  of  a  sphere  a  distance  h  above  the  ground 
plane.  The  surface  of  the  sphere  is  assumed  to  be  a  diffuse  isotropic 
radiator;  if  the  line  segment,  R,  connects  the  field  point  and  an  ele¬ 
ment,  dA,  of  the  surface  of  the  sphere,  then  the  radiation  from  dA  in 
the  direction  of  the  field  point  is  proportional  to  cosy,  where  y  is 
the  angle  between  R  and  the  outward  normal  to  dA.  The  radiation  from 
dA  which  is  intercepted  by  a  detector  placed  at  P  is  proportional  to 
cosB,  where  3  is  the  angle  between  R  and  the  normal  to  the  detector. 

It  is  required  that  ok  (tt/2);  for  o>(it/2)  the  detector  truncates  the 
sphere  to  less  than  a  hemisphere. 

Using  the  notation  of  Figure  12,  the  total  flux  incident  on  the 
detector  can  be  expressed  as: 


F(a,r,h,B) 


A 


cosy  dA 
~ rr  R7" 


(3) 


where 

-2  -1 

F  is  the  flux  at  the  detector  in  cal -cm  -sec 

_2 

S  is  the  spherical  surface  radiation  rate  cal-cm  -sec 
and  for  this  analysis  is  independent  of  the  position 
on  the  surface  of  the  sphere 

a  is  the  angle  between  r  and  the  detector  normal 


and 


B  is  the  radius  of  the  sphere  in  cm 

A  is  the  area  of  the  spherical  surface  which  is  above  the 
ground  plane  and  visible  from  P. 


To  evaluate  this  integral  the  various  line  segments,  angles,  and  inte¬ 
gration  limits  must  be  expressed  in  terms  of  the  constant  quantities 
a,r,h,B  and  the  integration  variables  0 , <j> .  These  expressions  include 
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Figure  12.  General  relationship  of  a  detector  at  P  to  a  truncated  sphere  at  C. 


R  =  (r2  +  B2  -  2rBcos©)^  (4) 

3  =  e(oc,e,4>)  (5) 

e  =  e(0)  (6) 

y  =  y(e,0) ,  y  =  0  +  e  (7) 


for  the  integration  limits 

0  <  ©  <  GM  (8) 

0M  =  cos-1  (B/r)  (9) 


0  <_  <_  4>q  (10) 

cb=T7if0<0  (11) 

o—o 

The  remaining  relations  are  developed  in  the  following  discussion. 

The  spherical  triangle  relation  can  be  used  to  express  3  in 
terms  of  other  variables  as 


cose  =  cosacose  +  sinasinecos<)>  .  (12) 

Forming  a  pair  of  right  triangles  from  the  triangle  which  contains  the 
angles  e,9,  it  is  seen  that 

Rsine  =  Bsin0  .  (13) 

The  small  diagram  in  the  upper  left  of  the  figure  is  used  to  obtain  one 
relation  for  the  line  segment,  y,  that  is: 

y  =  BsinOcos(iT-cj)0)  (14) 
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where  Bsine  is  the  radius  of  the  partial  ring  of  elemental  areas  at  e0<_ 
©<9^,  between  (-4>0 ,4>0 ) .  From  the  larger  diagram,  the  relation  for  y  is: 

y  =  tans(r  -  Bcos0)  (15) 


where 


6  =  sin  ^h/r). 

Equating  these  expressions  for  y  yields: 


(16) 


<f>0  =  cos  1 1  -tan6(r  -  Bcoso)/(Bsino)  j  ,  (17) 

which  expression  is  the  integration  limit  on  the  <f>  integral  for  0  <0s0. 


Using  the  triangle  formed  by  (P,C,G),  it  is  seen  that 


6  +  0+  cos  ^h/B)  =  tt/2 


or 


0Q  =  cos-1(h/r)  -  cos-1(h/B). 

The  expression  for  the  flux  at  the  detector  becomes 


=  f/Y 


*0  B2 

cosBcosy - sinOd^dO 


o  -<|> 


(18) 

09) 


(20) 


where 


♦o  = 


O<0<Go 


cos- 1|- tan  [sin-1  (h/r)]  (r  -  Bcos0)/(Bsin0)lo<p<p 


(21) 


e0  =  cos  1  (h/r)  -  cos  1 (h/B) 


(22) 
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(23) 


Gm  =  cos-1 (B/r ) 

v  =  0  +  e  (24) 
e  =  sin_1(Bsine/R)  (25) 
3  =  COS  1  (cosacose  +  sinasinecos<j>)  (26) 
R  =  (r2  +  B2  -  2rBcosQ)^  (27) 


2.2.2  The  Limit  of  the  Tangent  Sphere  (h  =  B) 

Before  proceeding  to  the  case  of  a  surface  burst,  it  is  worth¬ 
while  to  note  the  limiting  case  of  the  tangent  sphere  where  h  =  B; 
the  point  source  solution  should  result.  Setting  h  =  B  the  integral 
simplifies  to 


F  = 


_S 

TT 


Ej 

cosycosB^rr  sinedcfide. 
R 


0  -IT 


(28) 


Substituting  y  =  B/r  and  a,  -  cos"  y  and  performing  the  integration 
over  <f>,  the  intermediate  result  is 


■*/ 


COS  Y 

COSyCOSaCOSe 


2tty 


(1+Y  -2yCOSO) 


si  node. 


(29) 


Substitution  for  and  cosO  =  x  yields  the  algebraic  integral  form 

l 


F  =  2Sy  coset 


/- 


2  2 

[-y  +  (1  +  Y  )x  -  yX  ]  dx 

(1+Y  -  2yx) 


(30) 
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which  results  the  relation 

F  =  Sy2cosa  =  SB2  — ^  ,  (31) 

r 

which  is  the  point  source  solution,  as  expected. 

2.2.3  The  Limit  of  the  Hemisphere  (h  =  0). 

The  limiting  case  of  direct  interest  to  this  effort  is  that  of 
a  hemisphere  wherein  h  =  0;  in  which  case  the  general  integral  becomes 


F 


S 

IT 


COSyCOSB 


~  sin0d(j)d0 
R 


(32) 


Performing  the  ^-integration  and  separating  the  terms  with  the  factors 
cosa  and  sina,  the  integral  becomes 

0„, 


M 

/ct 


B 


F  =  Scosa  /  COSyCOSe  sin0d@ 

R 


a, 


+  S  —  sina  f 

’  / 


B  .  , 

cosysine  —  sin0d0  , 
R 


(33) 


which  may  be  expressed  as 

F  =  Scosal i  +  S  —  sinal2>  (34) 

where  Ix  and  I2  are  treated  separately  in  the  following 
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With  the  substitutions  y  =  B/r  and  x  =  cos0,  the  integral  Ix 


becomes 


[~y  +  (1  +  y^)x  -  yx2)  dx 

(1  +  y2  -  2yX ) 2 


or 


and  the  expression  for  F  becomes 
p  =  (SB2)  cosa  +  2S 


* 


I2sina 


(35) 


(36) 

(37) 

(37) 


The  substitutions  y  =  B/r  and  x  =  cos©  in  integral  I2  result 
in  the  form 


I2 


(x  -  xl  (1  - 

(1  +  y2  -  2yx)2 


dx. 


(33) 


2 

The  substitution  y  '=  1  +  y  -  2yx  is  helpful  in  evaluating  I2,  which 
becomes,  after  some  manipulation, 


I2  =  (sin_1y  -  y  (1  -  yY2).  (39) 

3 

Expansion  of  both  terms  demonstrates  that  I2  is  of  the  order  y  which 
is  anticipated  from  the  original  form  of  the  integral. 


*At  this  point  it  is  to  be  noted  that  for  a  =  0,  the  flux  at  the 
detector  is  one  half  that  of  the  tangent  sphere  solution,  which  is 
to  be  expected  from  the  symmetry  for  a  =  0. 
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Thus  the  flux  at  a  detector  at  field  point  P  for  a  radiating 
hemisphere  is  given  by 


r  SB  cosa  ,  S  .  )  .  -1  B 
F  =  -s—  — —  +  —  sina  <sin  — 

r  it 


which  for  B  «  r  is  approximately 


B 

r 


r 


,  (40) 


p  ~  SB  /  cosa  +  2sina  _B_ 


(41) 


2.2.4  Comparison  of  a  Special  Case  for  the  Hemisphere  with  Prior  Work. 

A  special  case  of  the  above  was  solved  earlier  by  a  number  of 
investigators  for  the  general  case  of  a  truncated  sphere  and  a  detector 
whose  normal  is  perpendicular  to  the  ground  plane.  For  this  case  the 
flux,  F',  at  the  detector  when  h  =  0  is  given  (reference  11  or  12)  in 
our  notation  as: 


F" 


.  -l 


sin 


'/(r/B)2  -  1 

(r/B) 


V(r/B)2~-~1) 

ir(r/B)2 


) 

I 


,  (42) 


Converting  this  form  to  one  similar  to  our  analytical  solution,  recog¬ 
nizing  a  =  tt/2,  yields  the  comparison: 


The  last  terms  are  seen  to  be  identical,  and  by  manipulating  the  first 
two  terms  in  the  relation  for  F'  it  is  seen  that  the  relations  for 
F  (a  =  tt/2)  and  F1  are  identities,  as  they  should  be. 
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2.2.5  Effect  of  Fireball  Radius  on  the  Geometrical  Shape  Factor. 

For  the  case  just  analyzed,  i.e.,  the  hemispherical  fireball 
radiating  to  a  detector  on  the  ground  whose  normal  is  perpendicular  to 
the  ground.  Figure  9  shows  the  strong  dependence  of  the  radiation*  on 
the  ratio  of  the  fireball  radius  to  the  detector's  range  from  the  burst 
point.  The  limit  3=0  corresponds  to  a  point  source  on  the  ground, 
which  will  not  give  any  radiation  perpendicular  to  the  ground  surface. 

The  radiation  arriving  at  the  same  detector  oriented  towards  the  fire¬ 
ball  does  not  depend  on  the  fireball  radius  at  all.  In  this  case  the 
flux  will  be  constant  no  matter  what  the  fireball  radius  is,  as  long 

as  the  power  radiated  from  the  whole  fireball  doesn't  change.  However, 

6  ^ 

the  geometrical  shape  factor  as  defined  will  be  equal  to  -y  .  It  is 

suggestive,  therefore,  to  also  present  the  dependence  of  another  quantity, 

1/ 

call  it  6,  where  8  =  —  .  Since  S  was  defined  as  the  surface  radiation 

32 

flux  from  the  hemispherical  fireball  surface,  the  power  radiated  out,  p, 
is  simply  2ttSB2.  Therefore  <5  is  proportional  to  the  flux  arriving  at  the 
detector  over  the  power  being  emitted,  whence  8  is  a  meaningful  quan¬ 
tity.  Figure  14  shows  the  dependence  of  8  on  3.  This  plot  shows  that 
when  the  fireball  is  nearly  on  top  of  the  detector,  the  orientation  of 
the  detector  between  these  two  limits  has  no  effect  on  the  received 
flux.  At  the  other  extreme,  the  smaller  the  fireball  radius  the  greater 
the  effect  of  detector  orientation. 

2.2.6  The  Fireball  Radius  Model 

Since  the  radiation  received  at  a  detector  on  the  ground  can 
be  a  strong  function  of  the  ratio  of  the  fireball  radius  to  detector 
ground  range  from  ground  zero,  a  reasonable  model  of  the  fireball 
radius  was  desirable.  Initially,  for  validating  the  computer  program 

*Shown  as  a  geometrical  factor,  that  is  equal  to  the  flux  (or  fluence) 
crossing  the  detector  divided  by  the  flux  (or  fluence)  emitted  from 
the  fireball's  surface. 


44 


1 


JL 


1 


10-1  - 

10"2  - 

/  DIRECTION  OF 

/  DETECTOR  NORMAL 

6 

/  POINTED  AT  - 

1  BURST  POINT 

10-3- 

POINTED  UP 

(I.E.,  NORMAL - 

TO  GROUND) 

10"4  - 

-5 

10  0  - 

0  0.2  0.4  0.6  0.8  1.0 

3 


Figure  14.  Variation  of  5  with  3. 
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code,  the  LAMB  fireball  model  was  used,  although  the  AFWL  never  intended 
it  to  represent  the  fireball  radius  accurately.  To  obtain  more  accurate 
thermal  radiation  estimates,  another  fireball  radius  model  was  required. 

Such  a  model  was  obtained  using  SPUTTER  data  available  on  the  FB21 
usertape.*  These  data  were  in  the  form  of  snapshots  at  certain  times. 

Since  the  early  time  fireball  was  difficult  to  identify  from  these  data, 
the  LAMB  fireball  model  was  used  for  the  early  time.  The  LAMB  fireball 
is  only  slightly  smaller  than  the  estimated  radius  up  until  12  milli¬ 
seconds  when  they  are  equal.  For  times  beyond  12  milliseconds  the  LAMB 
fireball  is  certainly  too  large;  and  the  model  developed  for  this  effort 
is  much  better.  The  new  representation  of  the  fireball  radius  should  be 
valid  up  until  about  600  milliseconds  for  a  1  kiloton  free-field  environ¬ 
ment.  It,  therefore,  provides  data  until  750  milliseconds  when  scaled 
to  the  1  kiloton  surface  burst. 

The  fireball  model  used  in  this  effort  is  presented  here,  and 
is  defined  for  scaled  time,  t,  where 

1  /3 

x  =  t/W  ,  W  =  yield  in  kilotons,  t  =  time  after  burst  in  seconds. 


When  x  <  0.012825 


R  =  2.568  x  102  t  0,395 


1 .6235  x  104t3,4962, 


when  0.012825  <  x  <  0.1 


R  =  1 .1565  x  102  x 


0.2119 


when  0.1  <  x  <_  0.2 

R  =  9.703  x  101  x  °'1357,  and 


*A1 though  only  three  Sputter  calculations  performed  at  the  AFWL  have 
been  documented  (references  4  through  6),  magnetic  tapes  were  distri¬ 
buted  throughout  the  community  that  contain  the  detailed  results  for 
each  calculation.  These  were  called  usertapes . 
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Table  2.  SPUTTER  data  from  FB21  usertape 
(1  kt  free-field  at  sea  level) 


ESTIMATE 
BOUNDS  FROM 

OF  FIREBALL  RADIUS 
SPUTTER  DATA  (FB21) 

FIREBALL 

RADIUS 

MODEL 

(m) 

TIME 

(MILLISECONDS) 

LOWER 

(m) 

UPPER 

(m) 

1 

19.2 

19.6 

16.8 

5 

35.0 

35.1 

31 .7 

10 

42.2 

42.4 

41 .6 

15 

46.8 

47.1 

47.5 

20 

50.3 

50.6 

50.5 

30 

55.0 

55.4 

55.0 

40 

58.8 

59.4 

58.5 

50 

61 .5 

62.2 

61  .3 

60 

64.0 

64.8 

63.7 

70 

65.3 

66.2 

65.8 

80 

67.4 

68.3 

67.7 

90 

68.8 

69.6 

69.4 

100 

71.2 

72.0 

71 .0 

150 

73.1 

75.4 

75.0 

200 

76.0 

79.2 

78.0 

300 

76.0 

79.5 

78.0 

400 

77.9 

77.9 

78.0 

500 

78.4 

78.4 

78.0 

600 

76.3 

76.3 

78.0 

700 

74.5 

74.5 

(not  modeled) 
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when  0.2  <  x  <  0.6 


R  =  78 

where  R  is  the  scaled  radius  of  the  fireball  in  meters/Kt^3.  The 
radius  of  the  fireball,  r,  is  then 

r  =  RW1/3 

1  /3 

For  scaled  times  beyond  0.6  seconds/Kt  this  fireball  radius  model 
is  not  defined. 

For  future  reference,  Table  2  presents  estimates  interpolated 
from  the  SPUTTER  usertape  data.  Where  the  fireball  edge  criteria 
clearly  defines  the  edge  the  bounds  will  be  equal.  The  values  for  the 
fireball  radii  were  taken  at  cell  boundaries,  and  in  some  cases  the 
differences  between  the  lower  and  upper  bound  on  the  radius  is  a  cell 
width.  Also  shown  for  comparison  are  the  results  of  using  the  model 
derived  from  these  data.  A  plot  of  the  fireball  radius  versus  time 
resulting  from  use  of  this  model  can  be  found  in  Section  4. 

2.3  RESULTS  FROM  SAI  THERML  CODE 

The  results  presented  in  this  section  as  well  as  those  in  the 
appendices  represent  estimates  of  nuclear  radiation  environments  from 
a  1  kiloton  surface  burst.  The  surface  burst  is  represented  as  a  1 
kiloton  hemisphere,  and  therefore,  the  radiation  output  has  been  com¬ 
puted  by  scaling  out  to  a  2  kiloton  free-field  (but  at  sea  level)  air 
burst.  No  account  has  been  made  of  any  contribution  due  to  scattering 
from  the  ground  or  other  objects.  The  ground  is  represented  as  a  per¬ 
fect  absorber  of  radiation  falling  upon  it,  except  for  the  radiation 
forming  the  fireball. 


49 


If  the  ground  is  in  fact  a  perfect  reflector  then  yield  doubling 
of  the  radiative  output  from  the  fireball  is  exact.  However,  when  the 
assumption  regarding  scattering  from  the  ground  is  invalid  it  can  lead 
to  significant  errors  for  some  cases.  If  on  the  other  hand,  the  earth 
is  in  fact  a  perfect  absorber  then  the  yield  doubling  is  inappropriate. 
The  appropriate  choice  for  airblast  is  to  double  the  yield.  The  approach 
taken  here  for  radiation  is  to  treat  the  ground  as  a  perfect  reflector 
inside  the  fireball,  but  then  assume  it  absorbs  all  radiation  incident 
upon  it  outside  of  the  fireball.  The  effect  of  these  assumptions  on 
the  calculation  of  flux  and  fluence  on  some  detectors  can  be  substantial. 
The  results  for  a  detector  whose  normal  is  elevated  90°,  i.e.,  the 
detector  is  lying  on  the  ground  facing  up,  will  obviously  not  be  affected 
by  the  assumption  about  how  the  earth  scatters.  However,  the  results  for 
such  a  detector  can  be  a  strong  function  of  the  fireball  radius  for  small 
radii  (see  Figure  13),  and  the  accuracy  of  the  fireball  radius  depends 
upon  the  appropriateness  of  yield  doubling.  On  the  other  hand,  the 
detector  whose  normal  is  at  elevation  0°  (facing  ground  zero)  will  not 
depend  directly  on  the  fireball  radius  but  will  be  influenced  by  the 
yield  doubling  models  in  terms  of  the  timing  and  magnitude  of  the  power 
out  as  well  as  scattering  from  the  ground.  The  effect  of  the  approach 
taken  will  be  discussed  further  in  Section  2.4. 

The  SAI  THERML  code  was  used  to  predict  the  environments  inci¬ 
dent  on  a  detector  at  ground  level.  The  ranges  considered  were  30.48, 
60.96,  91.44,  152.4,  304.8,  and  609.6  meters  from  ground  zero.  The 
results  shown  in  the  appendices  present  the  radiation  reaching  the 
detector  until  either  the  fireball  sweeps  over  the  detector  or  750  mil¬ 
liseconds  has  past.  Recall  that  this  latter  time  represents  the  limit 
(when  scaled)  on  the  model  for  the  fireball  radius.  For  all  results 

shown  the  fluence  is  given  in  cal /cm  ,  and  both  the  flux  and  average 

2 

flux  are  given  in  cal /cm  /sec. 
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Before  discussing  the  results  found  in  the  appendices,  it  is 
worthwhile  to  show  the  geometry  as  well  as  the  effect  of  the  orienta¬ 
tion  of  the  detector  on  some  of  the  results  for  an  interesting  ground 
range.  Figure  15  shows  the  detector  orientation  in  terms  of  an  eleva¬ 
tion  angle,  e,  measured  in  degrees.  The  normal  to  the  detector's 
surface  is  seen  to  be  in  the  plane  normal  to  the  earth  that  passes 
through  ground  zero  and  the  field  point.  Figure  16  traced  from  the 
appendices  shows  the  effect  of  detector  orientation  on  flux  from  a 
1  Kt  surface  burst  crossing  a  detector  152.4  meters  from  ground  zero. 
Figure  17  shows  the  effect  on  fluence.  Four  orientations  of  the 
detector  normal  were  considered:  0,  30,  60  and  90°  orientation. 

Results  for  30°  are  not  shown  since  they  are  almost  coincident  with  the 
0°  results.  As  expected,  for  detectors  facing  up,  i.e.,  90°  orientation, 
the  flux  and  fluence  environments  are  significantly  reduced.  Additional 
curves  similar  to  these  can  be  obtained  by  tracing  over  the  appropriate 
ones  in  the  appendices.  Appendix  A  contains  the  results  for  0°,  B  for 
30°,  C  for  60°  and  D  contains  the  90°  results. 

Some  of  the  results  are  summarized  in  Figures  18  through  21 
which  together  show  the  flux  and  fluence  at  each  ground  range  for  two 
detector  orientations  (0°  and  90°).  The  data  for  30.48,  60.96  and  91.44 
meters  ranges  include  times  up  until  engulfment  by  the  fireball.  Beyond 
these  times,  detailed  radiation  transport  solutions  are  required.  Al¬ 
though  these  data  are  available  from  the  SPUTTER  usertapes,  they  were 
not  needed  for  this  effort.  Since  the  fireball  model  does  not  predict 
growth  beyond  98.15  meters,  the  data  for  152.4,  304.8  and  609.6  meters 
are  complete  out  to  750  milliseconds. 

It  can  be  useful  for  purposes  of  simulation  to  have  estimates 
of  the  thermal  environment  at  the  time  of  shock  arrival.  Table  3  pre¬ 
sents  the  overpressure  at  time  of  arrival  based  on  a  model  by  Brode 
for  a  yield  of  2  kiloton  free-field  (or  a  1  Kt  surface  burst).  The 
equations  utilized  were 
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?7  and  ground  zero  in  same  plane 


Figure  15.  Detector  orientation. 
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Figure  16.  Effect  of  detector  orientation  on  flux. 
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Figure  17. 


Effect  of  detector  orientation  on  fluence. 
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Figure  18.  1  KT  surface  burst  :  detector  facing  GZ. 
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Figure  19.  1  KT  surface  burst  :  detector  facing  GZ. 
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Figure  20.  1  KT  surface  burst  :  detector  facing  up. 
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Figure  21.  1  KT  surface  burst  :  detector  facing  up. 
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Table  3.  Overpressures  and  arrival  times. 


GROUND 

RANGE 

(meters) 

FREE-FIELD 
OVERPRESSURE 
(psi ) 

SHOCK 

ARRIVAL 

TIME 

(ms) 

FIREBALL 

ARRIVAL 

TIME 

(ms) 

30.48 

3638 

2.5 

3.2 

60.96 

479 

14.3 

21. 

91  .44 

168 

37.0 

148. 

152.4 

49 

112. 

* 

304.8 

11 

409. 

* 

609.6 

3.7 

1190. 

* 

*Fireball  doesn't  grow  past  98  m. 
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(msec)  (45) 


t  (r,w) 

a 


0.5429w-21  ,185rw2/3+361 .8r2w] /3+2383r3 

w2/3+2.048rw1/3+2.687r2 


AP(t  ,w)  = 

a 


1 4843w 


1/3 


0.0135w1/3+t. 


i w2/3+0 . 671 5W1 /3t  +0 . 00481 3t  2 

_ 6. _ a 

w2/3+l .8836w1/3t  +0.0261 t  2 
a  a 


(psi)  (46) 


These  equations  were  obtained  from  reference  13,  page  25,  for  the  special 
case  t  =  tfl.  The  fireball  arrival  times  were  computed  with  our  model. 
These  arrival  times  were  then  used  with  the  data  in  the  appendices  to 
to  obtain  the  approximate  values  of  the  fluence  and  flux  at  the  time  of 
shock  arrival.  These  are  shown  as  Table  4.  The  data  do  not  exist  for 
the  609.6  meters  case  because  the  shock  arrives  at  a  time  later  than 
the  limit  on  the  fireball  radius  model. 


The  results  presented  in  the  appendices  provide  a  compendium 
that  can  be  used  to  interpolate  to  needed  conditions. 


2.4  UNCERTAINTIES  IN  THE  ESTIMATES  OF  NUCLEAR  THERMAL  RADIATION 

The  results  that  have  been  presented  so  far  are  estimates  of 
thermal  radiation  to  be  used  for  test  planning.  There  are  uncertainties 
associated  with  some  of  them.  Specifically,  there  exist  uncertainties 
associated  with  the  details  of  the  weapon-ground  and  thermal  radiation- 
ground  interactions.  These  uncertainties  have  not  been  quantified  as 
yet,  primarily  due  to  the  difficulty  of  such  a  task.  The  uncertainties 
introduced  by  the  modeling  of  SPUTTER  results  could  be  quantified;  how¬ 
ever,  it  is  beyond  the  scope  of  this  effort  to  do  so.  References  3 
through  5  provide  some  comparisons. 

The  uncertainties  that  are  discussed  here  should  be  borne  in 
mind  when  attempting  to  plan  what  environment  is  being  simulated. 
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Table  4.  Approximate  values  of  fluence  and  flux  at  shock  arrival. 


GROUND  1  KILOTON  SURFACE  BURST 

RANGE 

(meters)  FLUENCE  (CAL/ CM* )  FLUX  (CAL/Cfr/SEC) 


ORIENTATION  0°  30.48 

60.96 
91  .44 

152.4 

304.8 

ORIENTATION  30°  30.48 

60.96 
91  .44 

152.4 

304.8 

ORIENTATION  60°  30.48 

60.96 
91 .44 

152.4 

304.8 

ORIENTATION  90°  30.48 

60.96 
91  .44 

152.4 

304.8 


14 

1350 

10 

840 

35 

2100 

62 

360 

24 

11 

14 

1560 

9 

910 

37 

2230 

60 

360 

23 

10 

10 

1360 

7 

820 

28 

1780 

43 

270 

15 

7 

3 

790 

3 

500 

12 

840 

14 

100 

3 

2 
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2.4.1  Effect  of  Weapon/Ground  Interactions 

The  thermal  output  can  be  sensitive  to  details  of  weapon/ground 
interaction.  For  example,  the  thermal  radiation  can  be  shielded  by 
ejecta.  Further  there  may  exist  cooling  effects  on  the  fireball  due 
to  ground  mixing.  The  latter  changes  the  radiative  history  and  spectrum 
Reference  14  notes  a  dependence  of  the  radiative  history  on  height  of 
burst  for  several  NTS  shots.  For  example,  it  notes  that  there  appears 
to  be  strong  correlation  with  the  closeness  of  the  fireball  with  what 
amounts  to  be  the  thermal  partition.* 

SAI  modified  the  presentation  of  the  NTS  data  on  thermal  parti¬ 
tion  in  reference  14  by  plotting  it  versus  the  volume  of  a  sphere  that 
intersects  the  ground  and  represents  the  height  of  burst  and  fireball  ra 
dius  at  second  maximum  for  the  shots  indicated,  normalized  by  the  volume 
of  the  hemisphere.  This  yields  an  interesting  and  possibly  significant 
significant  trend  (Figure  22).  This  figure  suggests  that  within  the 
experimental  error  an  abrupt  change  in  the  character  of  the  thermal 
efficiency  occurs  when  the  weapon  is  very  close  to  the  ground.  Further, 
the  effect  seems  to  suggest  a  renormalization  with  a  radius  much  smaller 
than  the  fireball  at  second  maximum. 

SAI  immediately  though  of  the  x-ray  fireball.  The  SPUTTER 
radius  at  second  maximum  is  approximately  70  m  for  FB21 ,  1  kt  at  sea 
level,  whereas  the  early  time  fireball  is  only  4.5  m.  Assuming  this 
factor  is  appropriate  for  SMALL  BOY,  SUGAR  and  LACROSSE  just  for  the 
sake  of  argument,  SAI  obtained  a  figure  similar  to  Figure  23.  The 
implication  is  that  the  effect  is  pronounced  when  the  initial  x-ray 
deposition  (the  original  fireball)  is  striking  the  ground,  i.e.,  when 
ground  material  can  be  literally  blown  off  into  the  early  time  fire¬ 
ball  (and  hence  seriously  affect  its  development).  SAI  then  went  back 
to  the  SPUTTER  data,  estimated  the  x-ray  fireball  radius,  scaled  these 
for  the  nuclear  shots  shown,  and  then  used  these  to  normalize  the  data. 
The  result  is  presented  in  Figure  23. 

*The  apparent  thermal  radiated  yield  divided  by  the  total  yield. 
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Figure  22.  Thermal  efficiencies  (radius  at  2nd  max). 
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Figure  23.  Thermal  efficiencies  (radius  at  early  time). 
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Two  postulates  can  be  inferred  from  the  above.  First  if  the 
height  of  burst  is  much  greater  than  the  x-ray  burnout  region  the 
large  degradation  of  thermal  from  the  ground's  influence  (other  than 
ejecta  shielding)  will  be  relatively  small.  Secondly,  it  may  be  pos¬ 
sible  to  model  in  a  simple  way  from  existing  nuclear  data,  the  degrada¬ 
tion  of  thermal  radiation  efficiency  from  very  near  surface  bursts. 

The  data  in  reference  14  were  taken  from  cameras  stationed  on  the  ground, 
and  therefore,  the  correlation  may  only  hold  for  the  apparent  thermal 
efficiency  in  that  direction.  Since  this  is  the  direction  of  interest 
for  this  effort,  it  is  suggested  that  the  total  thermal  efficiency  com¬ 
puted  could  be  too  high  by  a  factor  of  7  for  an  actual  surface  burst. 
Although  the  correlation  may  quantify  the  effect  of  the  weapon-ground 
interaction  on  the  thermal  efficiency,  the  correlation  may  only  be 
fortuitous.  Further  investigation  is  required  and  believed  highly 
desirable. 

A  further  comment  is  appropriate  about  the  effect  of  weapon- 
ground  interactions  on  spectrum.  Reference  15  reports  that  near  sur¬ 
face  bursts  (as  well  as  below  surface  bursts)  modify  the  spectrum  of 
the  radiated  power  as  well  as  the  details  of  the  radiative  history. 

Finally,  the  effect  of  the  weapon-ground  interactions  on  fire¬ 
ball  development  has  not  been  quantified.  When  the  x-rays  penetrate 
the  ground  they  cause  vaporization;  the  vaporized  material  will  blow 
off  explosively  into  the  developing  fireball.  This  new  material 
will  affect  fireball  growth  and  its  radiative  history. 

In  summary,  the  degradation  of  the  radiative  history  from 
weapon-ground  interactions  has  not  been  quantified.  It  is  preferable 
that  the  curves  shown  in  the  appendices  be  used  to  define  the  simula¬ 
tion  allowing  study  of  synergistic  thermal  and  airblast  effects  for 
the  1  kiloton  surface  burst,  even  though  they  might  overpredict  the 
total  energy  radiated  out.  It  is  probable  that  such  a  simulation 
along  with  data  from  a  pure  airblast  environment  will  bound  the 
phenomena. 
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2.4.2  Effect  of  Thermal  Radiation/Ground  Interactions 

The  fireball  touching  the  ground  will  transfer  energy  to  the 
ground  by  a  combination  of  radiation,  conduction,  and  convection. 

Whether  these  interactions  significantly  affect  the  radiative  history 
of  the  fireball  is  still  an  unresolved  technical  issue.  It  is  not  anti¬ 
cipated  that  this  issue  will  be  resolved  for  some  time. 

2.4.3  Effects  of  Other  Uncertainties 

Crater  ejecter  can  obscure  the  fireball  somewhat  and  reduce  the 
amount  of  thermal  radiation  reaching  any  position  on  the  ground.  The 
effect  can  be  significant  for  low  yields.  What  effect  this  has  on  the 
radiative  output  has  not  been  quantified,  although  it  is  expected  that 
shielding  by  ejecta  will  not  be  as  important  as  the  effect  induced  by 
injection  of  ground  material  into  the  fireball. 
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SECTION  3 


RADIATIVE  OUTPUT  FROM  THERMAL  RADIATION  SIMULATION  DEVICES 

The  SAI  Thermal  Radiation  Simulation  (TRS)  devices  are  described 
briefly  in  this  section;  more  details  may  be  found  in  reference  16.  In 
addition,  this  section  presents  how  the  radiative  output  of  these  devices 
was  modeled  and  identifies  the  need  for  additional  characterization. 

3.1  THE  SAI  THERMAL  RADIATION  SIMULATOR 

During  1976,  a  preliminary  design  for  large  area  high  flux 
thermal  radiation  generator  was  conceived.  DNA  funded  a  small  program 
during  FY  77  to  test  this  concept.  The  results  of  the  initial  tests 
indicated  that  a  high  flux  (greater  than  200  cal/cm2-sec)  over  a  large 
area  could  be  achieved.  Subsequently  DNA  funded  a  full  scale  develop¬ 
ment  program  to  develop  a  TRS  device  which  would  produce  various  thermal 
radiation  pulses  of  high  flux  and  fluence  over  areas  exceeding  100  square 
meters.  It  was  to  be  self-consuming  upon  ignition  and  structurally 
simple. 

The  TRS  device  that  was  ultimately  used,  radiates  considerable 
amounts  of  energy  during  a  rapid  exothermic  oxidation  of  aluminum. 

Since  the  specific  energy  release  of  metallic  salt  is  relatively  high, 
most  of  the  chemical  energy  is  released  in  the  form  of  visible  and  near 
infrared  light. 

Though  many  pyrotechnic  mixtures  which  produce  intense  light 
are  available,  several  criteria  were  established  that  limited  the  choice 
of  mixture  and  led  SAI  to  adopt  aluminum  as  our  fuel.  These  criteria 
are : 

a.  The  mixtures  must  have  a  large  energy  release 
per  unit  weight  of  pyrotechnic  mixture. 

b.  The  pyrotechnic  mixture  must  be  safe  to  handle 
and  ship,  and  both  reactants  and  products  must 
be  non-toxic. 
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c.  The  mixture  must  have  a  low  cost  per  unit  weight. 

d.  The  material  must  be  readily  available  in  a  con¬ 
venient  form. 

Aluminum  powder  is  a  readily  available,  low  cost,  commercially  produced 
product  that  is  shipped  as  a  non-hazardous  substance,  and  upon  oxidation 
releases  high  energy  per  unit  weight. 

In  order  to  produce  high  intensity  bursts  of  thermal  energy, 
large  quantities  (>1.0Kg)  of  aluminum  must  be  ignited  and  burned  in  a 
relatively  short  time  (<1  sec).  To  accomplish  this  type  of  burn,  the 
metal  must  be  finely  divided  and  thoroughly  intermixed  with  oxygen. 

The  energy  release  associated  with  burning  aluminum  is  listed  in  Table  5. 


Table  5.  Aluminum  oxidation  reactions. 


Energy  Release 

Energy  Release 

per  gram  of 

per  Mole  of  Product 

Metal 

Reaction 

(cal/mole) 

(cal/gmM) 

2A1  +  -j-  02  -*  AI2O3 

3.89  x  105 

7.21  x  103 

Initial  work  conducted  during  1976  proved  that  aluminum  powder  could 
be  rapidly  mixed  with  oxygen  and  burned  to  created  high  flux  and  flu- 
ence  optical  pulses  over  large  areas.  Several  techniques  for  mixing  and 
igniting  the  A1  and  Op  mixtures  were  examined  during  the  tests.  The 
most  promising  technique  was  one  which  mixed  and  ignited  the  reactants 
by  use  of  a  small  quantity  of  high  explosive  (HE).  Figure  24  illus¬ 
trates  the  explosively  mixed  system  (E-system)  used  on  early  TRS  tests. 
Figure  25  shows  a  photograph  of  the  early  E-system  TRS  module. 
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Figure  24.  An  early  E-system  TRS  module 
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Figure  25.  Photograph  of  an  early  E-system  TRS  module. 


These  initial  tests  highlighted  areas  of  research  and  problems 
which  were  to  be  pursued  to  develop  a  useful  TRS.  Particular  problems 
and  research  areas  uncovered  in  the  early  TRS  work  were: 

a.  Early  TRS  systems  produced  excessive  airblast  during 
TRS  ignition. 

b.  Waveforms  of  early  TRS  systems  were  difficult  to 
control . 

c.  Scaling  laws  for  flux  and  fluence  generated  by  various 
yield  TRS  systems  had  to  be  developed. 

d.  Techniques  for  varying  pulse  signatures  had  to  be 
devel oped . 

e.  Techniques  for  the  control  and  calculation  of 
TRS  debris  had  to  be  developed. 


3.1.1  Laboratory  Experiments 

The  most  serious  problem  area  was  the  airblast  output  generated 
by  TRS  ignition  and  the  lack  of  control  over  TRS  radiation  pulse  shapes. 
A  small-scale  laboratory  research  program  was  conducted  to  study  these 
problem  areas  and  to  find  mechanisms  for  their  control.  These  experi¬ 
ments  indicated  that  the  macroscopic  combustion  rate  of  A1  was  respon¬ 
sible  for  generating  blast,  and  that  control  of  the  A1  combustion  rate 
could  lead  to  control  of  the  TRS  pulse  shape.  Variables  that  can  affect 
the  combustion  rate  are: 

a.  The  mixing  rate  of  A1  and  02  in  a  TRS, 

b.  The  A1  particle  size,  and 

c.  The  concentration  of  02  available  to  the 
Al-02  reaction. 
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Figure  26.  The  effect  of  the  O2/HE  weight  ratio  on  TRS 
fluence  of  E-system  TRS  modules. 
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A  research  program  was  conducted  to  examine  the  effect  of  these  as  well 
as  other  variables  on  the  thermal  radiation,  airblast,  and  debris  from 
the  simulator.  The  variables  considered  were: 

a.  Mixing  Technique.  High  explosive  dispersal 
(E-system)  and  high  pressure  gas  dispersal 
(C-system)  of  A1  into  an  oxygen  atmosphere 
were  examined. 

b.  Mixing  Rate.  The  mixing  rates  (i.e.,  the  high 

explosive  to  A1  weight  ratio  or  A1  weight  to 
driver  gas  pressure  ratio)  of  A1  and  CL  were 
varied.  ^ 

c.  Oo/Al  Ratio.  The  weight  ratios  of  A1  to  O2  in 
the  TRS  were  varied  to  examine  reactant  concen¬ 
tration  effects  on  TRS  output. 

d.  A1  Particle  Size.  The  size  of  A1  particles  and 
particle  size  distribution  in  a  mixture  of  A1 
powders  were  varied. 

e.  Physical  Size.  The  physical  size  of  systems  were 
varied  to  determine  A1  yield  scaling  effects  on 
TRS  output. 

The  approach  was  to  change  one  variable  while  holding  the  others 
constant.  Flux,  fluence,  and  airblast  were  measured  on  each  experiment. 
These  data  were  recorded,  analyzed,  and  used  to  design  larger  TRS  sys¬ 
tems.  Figure  26  represents  the  type  of  data  taken.  The  laboratory 
experiments  are  described  in  greater  detail  in  Reference  15. 


3.1.2  Small  Field  Experiments 

Since  TRS  systems  of  various  sizes  were  required  for  different 
thermal  radiation  experiments,  the  concept  of  modularized  TRS  units 
came  into  being.  To  create  larger  TRS  systems,  TRS  modules  are  added 
together.  Practical  dimensions  for  a  module  were  determined  by  study 
of  the  strength  of  different  inflated  plastic  cylinders.  TRS  modules 
with  length  to  diameter  ratios  of  5  and  inflated  to  an  overpressure 
of  0.1  psi  were  found  to  have  sufficient  structural  strength  to  resist 
deformation  by  35  mps  winds.  Further,  a  survey  of  common  plastic 
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materials  indicated  that  "lay-flat"  plastic  tubing  in  1  m  widths  (one 
meter  lay-flat  tubing  inflates  to  a  cylinder  with  a  diameter  of  0.64 
meters)  was  readily  available,  and  therefore  was  chosen  for  the  TRS 
module. 

Small  field  experiments  were  performed  to  create  a  controlled 
full-scale  TRS  module  and  to  examine  the  effects  of  A1  yield  scaling 
in  TRS  systems.  In  addition,  the  dimensions  of  these  tests  allowed 
examination  of  TRS  debris  problems  which  could  not  be  studied  on  smaller 
tests.  Field  experiments  were  conducted  at  the  explosive  test  facility 
leased  by  SAI  near  Pleasanton,  California.  A  series  of  12  tests  were 
conducted  during  this  portion  of  the  research  program.  TRS  flux,  fluence, 
and  airblast  output  were  measured  on  each  test.  Debris  and  TRS  fireball 
characteristics  of  each  experiment  were  recorded  with  high  speed  photo¬ 
graphy.  The  weight  ratios  of  A1  to  high  explosive  and  0^  to  A1  were 
held  constant  between  lab  and  field  experiments.  Aluminum  powder  yields 
on  these  field  scale  experiments  were  increased  by  a  factor  of  20  over 
the  lab  tests. 

3.1.3  Lab  and  Small  Field  TRS  Radiative  Output 

a)  Flux 

Figure  27  presents  the  time  dependent  power  output  of  lab  and 
field  TRS  systems.  Po wer  outputs  of  the  two  experiments  have  been 
scaled  to  equivalent  A1  powder  yields.  Variables  affecting  Al-0?  com¬ 
bustion  rates  have  been  held  constant  between  the  two  experiments. 

b)  Fluence 

Fluence  outputs  of  the  small  field  TRS  systems  followed  trends 
that  were  in  agreement  with  the  lab  experiments.  Field-scale  E-system 
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Figure  27a.  Lab-scale  E-system  TRS. 
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Figure  27b.  Small  field-scale  E-system  TRS. 


Figure  27.  Scaled  power  curves  for  two  E-System  TRS  designs 


TRS  designs,  for  given  C^/HE  weight  ratios,  produced  scaled  total  flu- 
ence  outputs  (cal/gm)  which  were  identical  to  those  produced  by  similar 
lab  experiments.  Field-scale  C-systems  TRS  designs  which  had  mixing 
rates  and  reactant  concentrations  similar  to  those  of  given  lab  systems 
also  produced  scaled  total  fluence  outputs  that  were  equivalent  to  lab 
experiments. 

3.1.4  Large  Field  Experiments 

The  purpose  of  the  large  field  experiments  was  to  examine  C-  and  E- 
system  multiple  module  effects  on  TRS  output.  Fourteen  large  tests 
were  completed  during  this  phase  of  the  research  program.  Five  events 
conducted  during  1976  were  also  used  as  data  sources.  With  the  comple¬ 
tion  of  the  large  field  experiments,  sufficient  data  had  been  compiled 
to  allow  completion  of  the  TRS  development  program. 

Four  different  experimental  set-ups  were  used  to  examine  large 
TRS  radaition  output.  The  set-ups  were: 

a.  Vertical  Single  Module  C-system  TRS  Designs.  A  vertical 
stand  which  held  a  single  A1  Fluidizer  unit  and  a  single, 
large  oxygen  bag  was  constructed  (not  shown). 

b.  Vertical  Multi-Module  C-system  TRS  Design.  A  vertical 
stand  which  held  4  A1  fluidizer  units  and  four  oxygen  bags 
was  constructed  (Figure  28). 

c.  Vertical  Multi-Module  E-system  TRS  Design.  A  vertical 
stand  which  held  12  E-system  TRS  modules  was  constructed 
(Figure  29) . 

d.  Horizontal  Multi-Module  E-system  TRS  Design.  A  horizontal 
frame  which  held  8  E-system  TRS  modules  and  could  be  raised 
to  a  height  of  6  meters  was  constructed  (not  shown). 

Set-ups  "a"  through  "c"  were  instrumented  with  FPT-100  photo¬ 
diodes;  CELESCO  LC  71,  PCB,  and  TYCO  HFG  airblast  gauges;  Medtherm 
GT-1600  high  speed  calorimeters,  SAI  fluence  calorimeters;  and  SAI 
aspirating  thermocouples.  Set-up  "d"  was  instrumented  with  the  CELESCO 
LC  71  airblast  gauge  and  a  Medtherm  GT-1600  highspeed  calorimeter. 
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Figure  28.  A  4  module  vertical  C-system  TRS  array. 
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Experiments  with  the  A1  fluidizer  and  a  cost  analysis  of  hard¬ 
ware  associated  with  it  indicated  that  the  fluidizer  should  inject  4  to 
5  kg  of  A1  powder  into  a  single  TRS  module.  Field  experiments  proved 
that  the  fluidizer  could  easily  spray  and  mix  A1  over  a  vertical  dis¬ 
tance  of  6  meters.  To  maintain  structural  integrity  in  the  module,  a 
plastic  cylinder  with  an  L/D  ratio  of  4  was  chosen  to  contain  the 
oxygen  of  the  TRS.  When  inflated  to  an  overpressure  of  0.1  psi,  the 
C-system  TRS  module  became  a  cylinder  with  a  length  of  6  meters  and  a 
diameter  of  1.52  m.  This  cylinder  provided  an  Al/O^  weight  ratio  in 
the  TRS  module  of  3.5.  Peak  power  outputs  for  the  single  C-system  TRS 
module  vary  from  1  x  107  cal /sec  to  3.5  x  107  cal /sec.  The  total  energy 
output  of  a  module  is  approximately  9  x  10^  calories. 

3.1.5  Large  Field  TRS  Radiative  Output 

a)  FI ux 

The  peak  power  generated  by  each  TRS  studied  in  this  portion 
of  the  research  program  was  measured  in  a  plane,  perpendicular  to  the 
plane  of  the  TRS.  In  general  ,  the  peak  power  generated  by  TRS  systems 
with  multiple  modules  of  similar  combustion  properties  (i.e.,  Al/Og  and 
Al/HE  weight  ratios,  mixing  techniques,  and  A1  yields),  was  a  linear 
function  of  the  number  of  modules  composing  the  TRS.  Figure  30  presents 
the  normalized  peak  power  outputs  of  several  multi -module  TRS  systems. 

The  peak  power  output  of  the  C-system  module  was  studied  on 
several  experiments.  In  general,  the  C-system  module  can  rise  to  peak 
power  over  controlled  time  periods  which  can  be  varied  from  25  msec  to 
500  msec.  The  peak  power  output  is  dependent  upon  the  TRS  ignition 
system,  fluidizer  mixing  rate,  and  A1  particle  size. 

b)  Fluence 

In  general,  two  problems  were  investigated  in  the  TRS  fluence 
study.  First  the  fluence  output  of  C-system  TRS  module  was  measured, 
and  second  the  effect  of  multiple  modules  on  TRS  fluence  output  was 
investigated. 
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Normalized  Peak  Power  (Power  of  Array  of  12/Power  of  1  Module) 
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Figure  30.  Normalized  power  outputs  of  multiple  module  TRS  arrays. 
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During  the  investigations  it  was  found  that  TRS  systems  which 
burned  18  y  A1  powder  were  more  efficient  than  TRS  systems  which  burned 
6  y  A1  powder.  For  given  combustion  conditions  (i.e.,  Al/C^  and  HE/C^ 
ratios,  ignition  conditions,  and  A1  yield)  the  scaled  fluence  output  of 
multi -module  TRS  systems  was  similar  to  the  single  module  output.  This 
indicated  that  the  energy  output  from  multi -moduled  TRS  systems  in¬ 
creases  linearly  with  the  number  of  modules  used  in  the  TRS.  However, 
more  needs  to  be  learned  about  the  transport  of  radiation  through  any 
intervening  TRS. 

With  the  completion  of  the  initial  TRS  testing  program,  a  large 
high  flux  and  fluence  TRS  had  been  developed.  Fluxes  of  200  cal/cm  -sec 
and  fluences  of  100  cal/cm2  had  been  achieved.  Since  the  completion 
of  the  research  program  the  TRS  has  been  used  to  successfully  irradiate 
soils,  aircraft,  shelter,  missiles,  and  tank  components.  Some  problems 
concerning  thermal  pulse  rise  time,  convective  vs  radiative  heating, 
and  pulse  reproducibility  have  been  encountered  since  the  initial 
research  project.  These  problems  need  to  be  further  addressed. 

3.2  COMPUTING  THE  RADIATIVE  OUTPUT  OF  TRS  DEVICES 

When  the  TRS  device  is  ignited  the  aluminum  near  the  source  of 
ignition  combusts.  The  chemical  reaction  spreads  in  a  complex  manner 
depending  on  where  the  aluminum  happens  to  be,  how  much  oxygen  is  avail¬ 
able  and  how  much  energy  is  available  to  initiate  the  combustion.  This 
process  quickly  forms  a  radiating  growing  fireball,  which  eventually 
cools.  Ideally,  the  combustion  occurs  homogeneously  and  uniformly 
throughout  the  fireball  simplifying  the  process  of  estimating  the 
radiative  output.  Practically,  this  doesn't  occur,  and  therefore, 
the  approach  taken  in  computing  the  radiative  output  is  approximate. 
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During  this  effort  a  computer  program  that  treats  each  module 
as  an  ensemble  of  point  sources  was  developed.  Each  module  is  treated 
as  a  linear  array  of  point  sources  along  the  module  centerline.  The 
program  has  been  written  to  facilitate  its  modification  to  treat  a 
module  as  a  surface  or  volume  radiator  as  appropriate.  If  treated  as  a 
volume  radiator,  some  geometrical  transmission  loss  factors  will  need  to 
be  incorporated. 

The  calculation  of  the  fluence  and  flux  are  straightforward 
when  the  thermal  radiation  devices  are  considered  as  an  ensemble  of 
point  sources.  The  fluence  F  is  simply 

cos a. 

F  =  £  S.  - J-  (47) 

i=l  4-rrr.j c 

where  S.  is  the  source  strength  of  the  i^  point  source  in  calories 

r.  is  the  vector  from  the  detector  to  the  i^ 
point  source 

a.  is  the  angle  between  r.  and  the  detector's 
1  1 

normal ,  n. 

N  is  the  total  number  of  point  sources  in  the 
ensembl e. 

If  one  device  is  represented  by  M  point  sources  then 

Sq 

S.  =  -jyj-  where  is  the  device's  total  output. 

Figure  31  shows  the  geometry  used  to  calculate  the  radiative 
output  from  the  point  source.  The  coordinate  system  chosen  here 
requires  the  xy  plane  to  be  the  ground,  and  the  xz  plane  should  con¬ 
tain  the  detector.  Then, 
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Figure  31 . 


Geometry  for  computing  contribution  from  k 


th 


source. 
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(48) 


cosa^  =  cosepcose^cose^,  +  sinepSine^ 


or 


cosepcose^  +  h^.sincp 


(49) 


whence 


^cosE^cose^  +  h^-sinsp 


} 


(50) 


The  representation  of  a  module  as  a  set  of  point  sources  (more 
than  one)  is  only  necessary  for  detectors  placed  close  to  it,  as  illus¬ 
trated  in  the  following  application  of  the  module  thermal  predictor 
program  to  the  configuration  in  Figure  32.  The  detector  is  oriented 
perpendicular  to  the  ground  (x-y  plane);  its  normal  along  x.  The  center- 
line  of  the  vertically-oriented  module  is  at  a  distance  R  from  the 
detector  at  an  angle  of  45°  (in  the  x-y  plane).  The  module  center  is 
placed  at  a  height  equal  to  one-half  its  length.  The  length  of  the 
chosen  module  is  6.7  m.  Since  its  diameter  is  approximately  1.5  m, 
ranges  closer  than  1  m  are  not  included.  The  module  was  represented  as 
a  set  of  N  points  equally  spaced  on  its  centerline.  The  results  are 
shown  in  Figure  33  for  four  cases:  N  =  1,  3,  5  and  7.  The  results 
are  insensitive  to  N  for  distance  greater  than  a  module  length  (6.7  m). 
For  distances  closer  than  this,  representing  the  module  as  1  point  is 
unacceptable.  A  good  rule  of  thumb  for  selection  of  the  number  of 
points  seems  to  be  to  set  N  equal  to  an  integer  obtained  by  normalizing 
the  length  of  the  module  by  the  range  of  interest  and  adding  1  to  the 
truncated  result. 


i .  e . ,  N 


(51) 
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Figure  32.  Orientation  of  detector  and  thermal  simulator  module. 
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Figure  33.  Variation  of  geometric  view  factor  with  number 
of  source  points  for  thermal  simulator  module. 
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Figure  34  presents  an  example  of  an  ensemble  of  eight  single 

thermal  radiation  simulator  modules.  For  this  ensemble  the  radiation 

can  be  computed  by  merely  summing  the  contributions  from  each  of  the 

single  modules.  An  example  of  such  computations  leads  to  Figure  35 

which  shows  the  effect  on  distance  away  from  the  ensemble.  At  1  m  the 

geometrical  view  factor*  is  0.041  per  square  meter.  The  output  of  one 

module  is  9  x  10^  cal,  whence  the  fluence  would  be  3.72  x  KT*  cal/m^ 

2 

or  37  cal /cm  .  The  points  shown  in  the  plan  view  are  each  in  the 
plane  which  is  orthogonal  to  all  the  module  centerlines  and  bisects 
each  of  them.  Figures  of  this  type  should  be  generated  for  reference 
purposes  once  the  thermal  radiation  simulator  output  is  better  char¬ 
acterized. 

3.3  NEED  FOR  ADDITIONAL  CHARACTERIZATION 

The  knowledge  of  the  radiative  output  of  the  thermal  radiation 
simulator  needs  to  be  improved.  One  aid  in  achieving  this  goal  is  to 
generate  empirical  data,  such  as  shown  in  Figure  36,  that  are  similar  to 
the  computational  data  shown  in  Figure  35.  Figure  36  shows  the  approxi¬ 
mate  maximum  view  factor  for  the  eight  module  array  shown  in  Figure  34. 
For  this  case,  iso-fluence  contours  are  obtained  by  multiplying  the 
view  factors  by  9  x  10  calories,  yielding  fluence  in  calories/rri  .  The 
results  of  Figures  35  and  36  will  be  readily  comparable  when  enough  is 
understood  about  the  physics  to  allow  simplified  models  to  be  generated. 
When  models  are  available,  the  fluence  at  various  detectors  and  orienta¬ 
tions  can  be  reliably  estimated.  Data  on  the  spectral  characteristics 
of  the  TRS  are  also  desirable  so  that  spectral  differences  between  nu¬ 
clear  environments  and  TRS  devices  can  be  taken  into  account.  At  this 
time  reliable  estimates  of  the  spectral  output  of  these  devices  are  not 
available. 


*Defi ned  by  setting  the  source  terms  of  each  module  to  unity  in 
equation  47.  This  statement  is  applicable  only  for  identical  modules. 
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R  =  distance  from  the  center 


of  a  given  module  to  the 


Figure  34.  Ensemble  of  thermal  radiation  simulator  modules. 
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Figure  35.  Effect  of  range  on  geometrical  view  factor. 
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Cross  Range  (m) 


Figure  36.  Contours  of  estimated  maximum  geometric  view  factors 
from  an  8-module  ensemble  (normalized  to  output  from 
one  module) . 
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SECTION  4 


TRS  PLACEMENT  IN  AN  AIRBLAST  ENVIRONMENT 

The  placement  of  the  thermal  simulation  devices  (TSD)  is  crucial 
to  an  adequate  simulation  of  the  thermal  radiation  environment.  Ideally 
the  easiest  placement  from  the  standpoint  of  producing  the  right  distri¬ 
bution  of  radiation  would  be  at  the  burst  location,  however,  since  the 
radiated  power  of  a  single  TSD  is  much  smaller  than  from  a  nuclear  device, 
a  large  number  of  them  are  needed  to  obtain  the  appropriate  levels.  For 
example,  assuming  superposition,  it  would  take  about  39,000  modules, 

g 

each  radiating  at  9  x  10  calories,  to  represent  a  1  kiloton  burst  having 
a  thermal  efficiency  of  35  percent.  A  typical  module  is  a  6  meter  high 
cylinder  of  1.5  meters  diameter.  Therefore,  if  these  modules  were  hexa- 
gonally  packed  they  would  fill  a  hemisphere  approximately  60  meters  in 
radius.  If  within  each  module  4  kilograms  of  aluminum  were  present,  over 
170  tons  of  aluminum  powder  would  be  consumed  within  this  hemisphere.  It 
is  worthwhile  noting  that  the  radius  of  a  one  kiloton  nuclear  fireball 
(surface  burst)  is  this  size  about  20  milliseconds  after  detonation. 

Figure  37  presents  the  relative  dimensions  of  the  initial  radius  of  a 
hemisphere  of  thermal  radiation  modules  required  to  produce  an  equivalent 
fluence  of  thermal  radiation  from  a  1  kiloton  of  the  stated  thermal 
efficiency  (35  or  5  percent)  with  appropriate  packing  (hexagonal  or 
square).  No  growth  data  are  shown  because  they  are  not  available. 

Figure  37  also  presents  the  time  of  arrival  from  Brode  (see  Section  2.3) 
and  the  fireball  model  for  1  kiloton  free-field  and  1  kiloton  surface 
burst.  It  should  be  noted  that  the  assumption  of  superposition  is  not 
believed  to  be  valid  in  these  circumstances.  It  is  expected  thalf the 
whole  ensemble  would  be  opaque  at  early  times  and  consequently  radiation 
will  be  issued  over  longer  periods  of  time,  probably  of  the  order  of  ten 
seconds.  This  duration  is  much  too  long  for  simulating  one  kiloton.  One 
can  expect  that  pulse  duration  will  scale  as  the  cube  root  of  the  ensemble 
yield. 
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Figure  37.  Fireball  radius-time  model  for  1  KT  surface  burst 
(and  1  KT  free-field). 
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The  feasibility  of  utilizing  such  a  large  ensemble  of  thermal 
radiation  simulators  has  not  been  investigated.  Considerable  work  would 
need  to  be  done  to  allow  their  use  in  such  an  arrangement.  Another 
approach  in  deploying  such  devices  is  to  substantially  reduce  their  dis¬ 
tance  from  the  field  point  of  interest.  The  thermal  radiation  can  then 
be  appropriately  modeled  at  the  field  point  of  interest,  although  it  will 
be  too  powerful  for  points  closer  to  the  device  and  too  weak  at  distances 
farther  away.  Careful  arrangement  of  an  ensemble  of  TSDs  can  also  be 
used  to  locally  improve  the  thermal  radiation  environment  simulation. 
Section  4.1  discusses  how  the  required  fluence  or  flux  might  be  obtained 
with  careful  arrangement  of  these  devices. 

Additional  effects  associated  with  use  of  these  devices  close  to 
targets  must  also  be  considered  when  arranging  these  devices  to  investi¬ 
gate  synergistic  blast  and  thermal  loads.  For  example,  care  must  be 
taken  not  to  appreciably  distort  the  blast  wave  flow  field  by  the  pres¬ 
ence  of  the  thermal  simulation  devices.  This  problem  as  well  as  tech¬ 
niques  for  its  solution  is  discussed  in  Section  4.2. 

Placement  of  these  devices  in  a  combined  experiment  also  leads 
to  consideration  of  the  effects  of  one  experiment  on  another.  This  is 
discussed  in  Section  4.3.  Finally,  additional  work  is  recommended 
before  fielding  these  devices  on  a  HE  test,  and  this  is  summarized  in 
Section  4.4. 

4.1  OBTAINING  THE  FLUENCE  AND  FLUX 

An  upper  bound  to  radiative  output  from  an  ensemble  of  thermal 
radiation  devices  can  be  estimated  by  multiplying  the  output  of  one 
device  by  the  number  of  devices.  If  each  device  is  considered  as  a 
point  source,  the  upper  bound  to  the  radiative  output  at  some  detector 
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can  be  estimated  by  placing  each  point  source  at  a  location  that  maxi- 

c  o  s  a 

mizes  the  geometric  view  factor,  -  g  ,  defined  earlier  in  Section  3. 

4  tit 

Similarily  the  maximum  possible  distance  an  ensemble  of  devices  can  be 

placed  to  obtain  some  fluence  can  be  estimated.  For  example,  if  the 

6  2 
output  of  one  device  is  9  x  10  calories,  a  fluence  of  10  cal /cm  is 

desired,  the  detector  is  oriented  to  obtain  the  maximum  fluence 

(i.e.,  cosa  =1),  and  ten  devices  are  to  be  used,  then  the  range  is 

simply 


r  = 


10  devices  x  9  x  10  calories/device 

(4tt)  10  cal/cm2 


846  cm 


8.5  meters 


Practically,  the  devices  will  not  all  be  positioned  at  the  same  loca¬ 
tion,  and  therefore  the  view  factor  for  some  of  the  devices  will  be 
less  than  its  maximum  possible  value.  Hence,  the  ensemble  will  need 
to  be  placed  closer  to  the  point  of  interest  and  the  estimated  radius 
would  be  further  refined  by  performing  a  set  of  calculations  (Section  3) 
and  interpolating  on  the  range.  One  effect  that  could  reduce  the  range 
further  would  be  obscuration  by  other  objects,  including  other  devices. 

Figure  38  is  used  to  illustrate  the  computing  of  the  total 

fluence  at  some  object.  Here  the  module  height  is  taken  to  be  6  meters, 

and  its  diameter  is  1.5  m.  In  this  example,  the  experimental  require- 
2 

ment  is  1  cal/cm  normal  to  the  object  3  m  above  the  ground.  Further 
it  is  desired  to  provide  a  window  of  width  W  through  which  the  airblast 
can  propagate  relatively  free  of  the  effects  of  the  TSD  fireball.  To 
simplify  the  problem  somewhat,  symmetry  is  employed  about  the  plane  con¬ 
taining  the  line  of  sight  and  normal  to  the  ground.  The  range,  R,  and 
angles,  9  and  <}>  at  which  each  ensemble  of  modules  is  placed  can  be  varied 
as  appropriate. 
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Figure  38.  An  example  of  module  placement. 
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Estimating  the  distance  the  modules  can  be  placed  from  the 


object 


6  -  2 
16  devices  x  9  x  10  calories/device  x  cose 

(4tt)  1  cal/cm2 


or  R  ~  34  meters  [cose] 


where  the  effect  of  different  <j>  has  initially  been  ignored.  If  the 
experimenter  were  to  orient  the  ensemble  such  that  £  was  90°  (see 
Figure  36,  which  also  showed  the  effect  of  various  orientations  on 
maximum  possible  fluence)*  then 


W  n  .  B 
~2~  =  Rsme  -  -y-  cose 

and  therefore  if  e  were  45°  R  would  be  24  meters  and  the  window 
width  would  be  25.4  meters.  If  instead  <f>  were  set  to  e  then  W  would 
be  increased  to  2Rs i n or  34  meters  for  the  same  e  of  45°. 

Estimating  the  output  of  those  same  two  arrangements  more 

2 

accurately,  the  fluence  is  found  to  be  1.36  cal/cm  when  the  center 

of  each  ensemble  is  a  range  of  24  meters  for  <f>  =  90°,  and 
2 

1.42  cal/cm  at  the  same  range  but  for  <j>  =  e.  The  comparison  shows 
that  the  experimenter  gains  by  using  the  latter  arrangement  in  addi¬ 
tion  to  providing  the  maximum  window.  The  only  loss  will  be  the 
possible  effect  of  transmission  through  adjacent  module  fireballs. 
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4.2 


SHOCK  INTERACTIONS  WITH  FIREBALL 


It  is  desirable  to  avoid  placing  the  thermal  radiation  simulators 
between  the  explosion  source  and  the  object  being  radiated.  The  result 
of  firing  the  thermal  radiation  simulator  module  is  to  create  a  hot  fire¬ 
ball  of  underdense  gas  that  is  eventually  in  pressure  equilibrium  with 
the  atmosphere  and  is  filled  with  the  products  of  combustion  (Al^O^, 
ignitor  residue,  and  small  pieces  of  mylar).  These  latter  constituents 
can  be  hot  particulates.  In  Figure  39,  the  propagation  of  a  shock  wave 
through  a  fireball  is  shown  qualitatively.  Since  the  fireball  is  hot, 
the  shock  will  propagate  through  it  more  quickly  than  through  the  sur¬ 
rounding  atmosphere.  It  will  also  compress  the  hot  fireball  and  move 
particulate  matter  within  it  closer  together.  This  will  tend  to  spread 
the  shock  out  as  shown  in  Figure  39,  introducing  two  dimensional  effects.* 

The  passage  of  a  planar  shock  through  a  hot  region,**  which  is 
in  pressure  equilibrium  with  its  surroundings  before  shock  arrival, 
appears  to  have  little  effect  at  high  overpressure  on  the  transmitted 
pulse.  This  has  been  demonstrated  using  two  calculations  performed  with 
the  SAI  one-dimensional  version  of  the  HULL  code  (References  17,  18). 
First,  a  reference  calculation  was  performed  at  sea  level  atmospheric 
pressure,  density,  and  internal  energy.  At  an  inlet  boundary,  a  nuclear 
airblast  waveform  was  used  which  corresponded  to  a  59.4  psi  peak  over¬ 
pressure  from  a  1  kiloton  surface  burst.  This  airblast  waveform  was 
obtained  using  the  AFWL  LAMB  code  (described  in  references  19  and  20). 

A  second  calculation  was  then  performed,  but  now  the  propagation  path 
contained  a  hot  region  starting  20  m  from  the  inlet  boundary  and  lasting 
for  3  m  as  shown  in  Figure  40.  Several  observer  locations  were  included 
in  each  calculation  in  order  to  allow  for  comparisons.  A  reflective 
boundary  was  placed  at  the  end  of  the  computational  mesh  so  that  signals 


*Since  the  fireball  is  a  finite  cylinder,  end  effects  will  lead 
quickly  to  three-dimensional  effects. 

**without  particulates 
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PLAN  VIEW 


SHOCK 


BEFORE  SHOCK  ARRIVAL 


Figure  39.  Airblast  interaction  with  TRS  device  fireball. 
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CLEAN  AIR 


INLET  CONDITIONS:  59.4  PSI  SHOCK  FROM  A  1  KT 

SURFACE  BURST  (USING  LAMB 
MODEL  FOR  WAVEFORM) 


FIREBALL  PRESENT 


REGION  IN  FRONT  OF  FIREBALL 

REGION  BEYOND  FIREBALL 

REFLECTIVE 

INLET— # 

*  * 

BOUNDARY 

BOUNDARY 

1 

Figure  40.  Idealization  in  one-dimension  of  shock  propagating 
through  fireball. 
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propagating  back  into  the  mesh  would  be  obvious.  The  spatial  ordinate 
is  placed  at  the  left  hand  side  of  the  figure,  140  m  from  the  1  kiloton 
surface  burst.  The  mesh  is  approximately  42  m  long,  therefore  placing 
the  right  hand  boundary  182  m  from  the  burst.  Both  calculations  were 
one-dimensional  rectangular  and  therefore  underpredict  geometrical 
attenuation.  However,  our  desire  here  is  to  assess  the  importance  of 
the  TRS  fireball  on  airblast  propagation.  Initial  conditions  for  each 
case  are  given  in  Table  6. 


Table  6.  Computational  parameters. 


Calculations  with  Fireball  Present 

P=P  ,  £  =  £  ,  1=1 ,1 00 

o  o 

p=pO/10,  e=10e  ,  1=101,115 

o 

p=p  ,  £=£  ,  1  =  116,215 

o  o 


Clean  Air  Calculation 
P=P  ,  £=£  ,  1=1 ,21 5 

O  0 

g 

where  £^=2. 08x10  ergs/gm 
P  =1 . 225x1 0"3gm/cm3 


IMAX  =215 
DX  =  20  cm 
XQ  =  1 .4  x  104  cm. 

OP  =  59.4  psi 

W  =  2  kilotons  (i.e.,  1  kiloton 
surface  burst) 

Tstart  '  93'65  msec 
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Figures  41  through  45  present  comparisons  of  the  two  calcula¬ 
tions  for  various  observer  locations.  The  observer  locations  are  shown 
in  Figure  40  at  5  m  forward  of  the  hot  fireball  and  2  m  and  7  m  beyond 
the  hot  region.  Each  calculation  began  94  milliseconds  after  detonation 
of  a  1  kiloton  device,  when  the  airblast  arrives  at  a  range  of  140  meters. 
The  wave  propagates  through  the  mesh  in  about  60  milliseconds.  The 
reflected  wave  does  not  arrive  at  any  observer  stations  shown  in  Figures 
41  through  45. 

When  the  shock  passes  through  the  hot  region  the  density  in  the 
fireball  is  increased  by  a  factor  of  about  15  to  20,  and  the  hot  com¬ 
pressed  region  travels  in  the  direction  of  the  incident  shock  wave.  The 
presence  of  this  fireball  causes  relatively  small  changes  to  the  trans¬ 
mitted  velocity  and  overpressure  peaks.  Specifically,  a  7.7  percent 
reduction  in  the  overpressure  peak  is  seen  2  meters  beyond  the  fireball 
as  well  as  a  5  percent  reduction  in  peak  velocity.  The  comparison  of 
velocities  at  this  observer  position  is  shown  as  Figure  41.  Further 
downstream,  7  meters  beyond  the  fireball,  the  peak  overpressure  is 
reduced  by  only  2  percent  and  the  velocity  by  about  one  percent,  with 
the  differences  diminishing  more  with  increasing  distance.  The  hot 
region  thus  causes  the  pulses  to  arrive  sooner,  together  with  a  broad¬ 
ening  of  the  pressure  and  velocity  peaks  especially  for  observers  close 
to  the  fireball.  The  fireball  also  modifies  the  upstream  results  as 
shown  in  Figures  44  and  45,  where  it  can  be  seen  that  overpressure 
drops  more  rapidly  and  velocity  more  slowly  at  late  times  with  the 
fireball  present.  These  calculations  do  not  include  the  effect  of 
any  hot  particulate  matter.  Since  the  ratio  of  total  particle  mass  to 
fireball  volume  is  approximately  the  same  as  the  ambient  air  density, 
the  effect  may  be  important  and  may  further  reduce  the  strength  of 
transmitted  pressure  pulse.  In  addition,  the  two-dimensional  effects 
have  not  been  estimated.  A  typical  calculation  that  might  be  done  to 
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Velocity  (cm/sec) 


I 


Figure  41 . 


Time  (sec) 


Comparison  of  velocity  at  observer  approximately 
2  meters  beyond  fireball. 
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Velocity  (cm/sec) 


Time  (sec) 


Figure  42.  Comparison  of  velocity  at  observer  approximately 
7  meters  beyond  fireball. 
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Overpressure  (psi) 


Time  (sec) 


Figure  43.  Comparison  of  overpressure  at  observer  approximately 
7  meters  beyond  fireball. 
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Overpressure  (psi) 


Time  (sec) 


Figure  44.  Comparison  of  overpressure  at  observer  approximately 
5  meters  in  front  of  fireball. 
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Velocity  (cm/sec) 


Time  (sec) 


Figure  45.  Comparison  of  velocity  at  observer  approximately  5  meters 
in  front  of  fireball. 


106 


estimate  the  importance  of  the  two-dimensional  effects  is  shown  as 
Figure  46.  An  inlet  boundary  would  again  be  used  to  introduce  an  air- 
blast  into  hydrodynamic  mesh,  and  a  reflective  boundary  would  be  used 
to  employ  symmetry.  Good  transmissive  boundaries  would  keep  the  costs 
down  by  not  including  more  mesh  than  needed. 

4.3  EFFECT  OF  OTHER  EXPERIMENTS 

When  designing  any  program  involving  more  than  one  experiment, 
the  interaction  of  experiments  must  be  considered.  For  example,  as  dis¬ 
cussed  previously,  the  presence  of  the  thermal  simulation  devices  may 
modify  the  airblast  at  the  target  of  interest.  Care  must  also  be  taken 
to  recognize  that  the  presence  of  an  ensemble  of  TRS  modules  can  also 
affect  the  airblast  incident  on  ancilliary  experiments.  The  magnitude 
of  these  effects  must  be  addressed  quantitatively. 

To  illustrate  these  points,  Figure  47  shows  a  possible  arrange¬ 
ment  of  three  experiments  being  performed  on  objects  A,  B,  and  C  using 
a  common  source  of  airblast.  The  angular  position  of  objects  A  and  B 
are  fixed  by  other  constraints,  and  their  ranges  are  set  by  required 
overpressure  levels.  C  is  assumed  to  be  in  the  sector  shown  bounded 
by  the  radials  to  A  and  B.  It  would  of  course  not  be  placed  to  close 
to  either  e  =  0°  or  3  =  a  so  as  to  effect  overpressures  at  A  or  B.  The 
radiative  environment  added  by  the  TRS  ensembles  associated  with  object 
C  can  be  calculated  at  positions  A  and  B  as  a  function  of  angle  3  to 
determine  the  best  location  for  C.  Similarly  the  environments  from  A 
and  B  must  be  considered  when  designing  experiment  C.  Shielding  on  the 
back  side  of  the  object  may  also  be  desirable  in  some  cases. 

In  summary,  the  TRS  modules  on  an  HE  event  will  involve  the 
introduction  of  another  variable  on  the  fidelity  of  the  simulation  on 
each  experiment. 
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INLET  BOUNDARY 


TRANSMISSIVE  BOUNDARY 


Figure  46.  An  example  calculation. 
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Range  (m) 


Figure  47.  An  example  experimental  layout. 
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4.4 


NEED  FOR  ADDITIONAL  WORK 


An  earlier  section  discussed  the  need  for  additional  work  to 
calibrate  the  radiative  output  of  a  thermal  radiation  simulator  module. 
This  section  summarizes  the  need  for  additional  work  so  that  these 
modules  can  be  used  in  an  ensemble  to  simulate  thermal  radiation  in 
airblast  environments.  In  this  summary,  it  is  assumed  that  the  exper¬ 
imenter  will  attempt  to  simulate  the  thermal  radiation  with  a  small 
ensemble  of  well  placed  modules  rather  than  using  a  large  hemisphere 
of  modules  near  the  HE  airblast  source.  Furthermore,  the  discussion 
is  restricted  to  ensembles  of  modules  where  each  modules  is  well  char¬ 
acterized. 

In  Section  3,  data  were  presented  (Figure  30)  that  indicated 
that  the  peak  power  radiated  was  only  a  function  of  the  number  of 
modules  in  an  ensemble  (data  limited  to  12  modules).  The  presence  of 
other  modules  does  not  seem  to  modify  the  pulse  length  and  radiative 
waveform  for  conditions  tested  so  far.  Unfortunately  these  comparisons 
are  not  well  documented.  In  addition  better  definition  of  the  trans¬ 
mission  losses  through  adjoining  modules  would  be  useful  and  would 
facilitate  modeling  of  the  thermal  environment.  This  work  needs  to  be 
done  in  order  to  determine  fluence  and  flux,  but  also  spectral  content. 
Whether  any  module  will  preferentially  absorb  radiated  energy  from  an 
ajoining  module  needs  to  be  known. 

The  dynamic  extent  of  the  fireball  must  also  be  known.  It  is 
known  that  the  burning  of  the  aluminum  powder  in  the  oxygen  environment 
leads  rapidly  to  a  fireball  that  does  grow.  Little  quantitative  data 
have  been  obtained.  The  effect  of  the  presence  of  other  modules  on  a 
fireball's  growth  also  needs  to  be  understood.  It  is  possible  that  a 
module  constrained  by  two  adjacent  modules  might  grow  faster  and  maybe 
further  in  the  direction  where  expansion  is  not  confined.  Figure  48 
presents  this  idea  pictorial ly  showing  a  plane  cutting  the  modules  per¬ 
pendicular  to  their  cylindrical  lines  of  symmetry. 
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Figure  48.  Possible  effect  of  confinement  on  fireball  growth. 
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Finally,  more  theoretical  and  experimental  work  are  needed  to 
better  define  the  effect  of  the  presence  of  the  remaining  TRS  fireball 
on  airblast  propagation. 
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SECTION  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


Although  techniques  are  now  available  for  rapidly  burning  alu¬ 
minum  in  an  oxygen  environment  in  order  to  produce  considerable  radiant 
energy  over  large  areas,  additional  data  are  needed  to  allow  proper  use 
of  these  techniques  when  simulating  the  thermal  radiation  environment  of 
a  nuclear  device.  Furthermore,  at  this  time  the  nuclear  environment  is 
only  partially  simulated  and  questions  remain  about  the  degree  of  simu¬ 
lation  fidelity. 

The  characterization  of  the  radiative  output  of  the  SAI  thermo¬ 
chemical  thermal  simulation  device  module  is  incomplete  in  a  number  of 
ways.  Although  some  theoretical  modeling  has  been  done  by  SAI  on  pilot 
small  scale  systems  (reference  21),  more  work  needs  to  be  done  for  recent 
large  scale  systems.  Data  that  are  essential  to  modeling  radiative  out¬ 
put  of  a  TDS  module  include: 

1)  The  fireball  extent.  After  TSD  module  ignition, 
the  TSD  fireball  begins  to  expand.  Data  on  its 
size  and  shape  from  initiation  to  stabilization 
are  needed  in  order  to  better  characterize  the 
TDS  as  a  source  of  radiation.  Some  data  exist  but 
need  to  be  analyzed  and  documented. 

2)  The  spectrum  of  radiation.  It  is  desirable  to 
obtain  spatial,  temporal,  and  spectral  data  on 
specific  intensity  from  the  fireball.  Some  inte¬ 
grated  data  exist;  however,  they  are  insufficient 
for  modeling  the  radiation  on  a  target  since  they 
are  not  resolved  in  frequency  and  space. 

Some  recent  experiments  that  have  been  performed  have  shown  that 
placing  one  TDS  module  behind  another  can  increase  the  radiation 
reaching  a  target.  The  increase  is  less  than  the  factor  of  two 
obtained  using  superposition.  The  effect  may  be  due  to  the  fact  that 
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the  fireball  nearest  to  target  is  partially  transparent;  however,  esti¬ 
mates  made  in  reference  21  imply  that  the  expected  fireball  is  much  too 
opaque  for  this  to  occur.  It  is  also  possible  that  the  adjacent  module 
constrains  the  other  fireball  from  radiating  in  its  direction  but  beyond 
its  position,  leading  ultimately  to  more  radiation  towards  the  target. 

This  needs  to  be  investigated  further  so  that  models  can  be  either  val¬ 
idated  or  improved  sufficiently  to  characterize  the  use  of  such  modules 
in  an  ensemble. 

There  exists  no  substantive  data  on  the  effect  of  the  TRS  devices 
on  propagation  of  air  blast,  although  very  recent  and  preliminary  air 
blast  measurements  on  Misers  Bluff  seem  to  indicate  that  the  air  blast 
propagating  through  a  TRS  attenuates  somewhat  as  the  one-dimensional 
model  indicates.  We  conclude  that  the  one-dimensional  effect  needs  to 
be  further  investigated  to  account  for  the  influence  of  hot  particles. 

The  importance  of  two-dimensional  effects  also  needs  to  be  analyzed. 

The  prediction  of  the  thermal  radiation  from  nuclear  bursts 
during  an  actual  war  requires  consideration  of  surrounding  conditions 
(such  as  clouds,  haze,  hills,  etc.),  as  discussed  in  reference  2.  Any 
description  of  the  environment  that  is  to  be  simulated  is  in  itself 
non-trivial.  Nevertheless  if  the  thermal  radiation  environment  from  the 
fireball  is  known,  these  other  effects  can  be  estimated.  Although  thermal 
environments  for  detonations  far  above  ground  can  be  well  characterized, 
the  use  of  free-field  thermal  radiation  descriptions  for  near-surface 
detonations  may  be  incorrect  both  in  flux  level  and  spectral  content. 
Radiation  out  to  space  can  be  much  different  from  that  reaching  the 
ground  for  near  detonations.  This  study  has  led  to  the  development  of 
a  correlation  between  thermal  efficiency  and  nearness  of  the  device  to 
the  ground.  Although  this  correlation  may  be  fortuitous,  it  warrants 
considerable  additional  investigation. 
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Finally  the  following  recommendations  are  offered  based  on 
efforts  conducted  under  this  contract.  In  the  area  of  thermal  radia¬ 
tion  simulation,  work  should  continue  in  order  to  develop  a  means  to 
achieve  repeatable  light-flux  variable-pulse-length  thermal  radiation 
waveforms  over  large  areas.  It  is  further  recommended  that  the  thermal 
radiation  environment  for  these  thermal  waveforms  be  well  characterized 
and  modeled.  In  addition,  scaling  laws  for  radiative  output,  full 
requirements,  and  number  and  placement  of  modules  should  be  determined 
for  the  simulators. 
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APPENDIX  A 


ESTIMATE  OF  THERMAL  RADIATION 
FROM  A  1  KILOTON  NUCLEAR  SURFACE  BURST 
DETECTOR  ELEVATION  0  DEGREES 


This  appendix  contains  plots  for  one  elevation  and  consists  of 
six  sets  of  three  figures.  Each  set  corresponds  to  a  particular  range: 
30.48,  60.96,  91.44,  152.4,  304.8,  or  609.6  meters  from  ground  zero. 

For  each  range  fluence,  flux,  and  average  flux  are  presented  as  a  func¬ 
tion  of  time.  Fluence  is  in  units  of  cal -cm-2,  flux  and  average  flux 
are  in  units  of  cal-cm-2-sec-1 . 


119 


FLUENCE  (CAL/CM  ) 


Figure  49.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  30.48  meters  ground 
range. 
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Figure  50.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  0  degrees  at  30.48  meters  ground 
range. 


121 


AVERAGE  FLUX  (CAL/Ci'i  /SEC) 


TIME  SEC 


Figure  51.  Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  30.48  meters  ground  range. 
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FLUEl\iCE(CAL/Chf) 


Figure  52.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  60.96  meters  ground 
range. 
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Figure  53.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  0  degrees  at  60.96  meters  ground 
range. 
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Figure  54.  Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  60.96  meters  ground  range. 
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Figure  55.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  91.44  meters  ground 
range. 
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FLUX  (CAL/ CM  /SEC) 


Figure  56.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  0  degrees  at  91.44  meters  ground 
range. 
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Figure  57.  Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  91.44  meters  ground  range. 
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Figure  58.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  152.4  meters  ground 
range. 
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Figure  59.  Flux  from  1  kt  surface  burst  crossing  detector 

wi~h  orientation  0  degrees  at  152.4  meters  ground 
range. 
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Figure  60.  Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  152.4  meters  ground  range. 
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Figure  61.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  304.8  meters  ground 
range. 
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Figure  62.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  0  degrees  at  304.8  meters  ground 
range. 
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Figure  63.  Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  304.8  meters  ground  range. 
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Figure  64.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  609.6  meters  ground 
range. 
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Figure  65.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  0  degrees  at  609.6  meters  ground 
range. 
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Figure  66.  Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  0  degrees  at  609.6  meters  ground  range. 
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APPENDIX  B 


ESTIMATE  OF  THERMAL  RADIATION 
FROM  A  1  KILOTON  NUCLEAR  SURFACE  BURST 
DETECTOR  ELEVATION  30  DEGREES 


This  appendix  contains  plots  for  one  elevation  and  consists  of 
six  sets  of  three  figures.  Each  set  corresponds  to  a  particular  range: 
30.48,  60.96,  91.44,  152.4,  304.8,  or  609.6  meters  from  ground  zero. 

For  each  range  fluence,  flux,  and  average  flux  are  presented  as  a  func- 
tion  of  time.  Fluence  is  in  units  of  cal -cm-  ,  flux  and  average  flux 
are  in  units  of  cal -cm-^ -sec-1 . 
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Figure  67.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  30  degrees  at  30.48  meters  ground 
range. 
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Figure  68.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  30.48  meters  ground 
range. 
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Figure  69.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  30.48  meters  ground  range. 
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Figure  70.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  30  degrees  at  60.96  meters  ground 
range. 
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Figure  71.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  60.96  meters  ground 
range. 
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Figure  72.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  60.96  meters  ground  range. 
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Figure  73.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  30  degrees  at  91.44  meters  ground 
range. 
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Figure  74.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  91.44  meters  ground 
range. 
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Figure  75.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  91.44  meters  ground  range. 
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Figure  76.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  30  degrees  at  152.4  meters  ground 
range. 
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Fiqure  77.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  152.4  meters  ground 
range. 
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Figure  78. 


Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  30  degrees  at  152.4  meters  ground  range. 
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Figure  79.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  30  degrees  at  304.8  meters  ground 
range. 
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Figure  80.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  304.8  meters  ground 
range. 
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Figure  81.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  304.8  meters  ground  range. 
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Figure  82.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  30  degrees  at  6G9.6  meters  ground 
range. 
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Figure  83.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  609.6  meters  ground 
range. 
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Figure  84.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  30  degrees  at  609.6  meters  ground  range. 
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APPENDIX  C 


ESTIMATE  OF  THERMAL  RADIATION 
FROM  A  1  KILOTON  NUCLEAR  SURFACE  BURST 
DETECTOR  ELEVATION  60  DEGREES 


This  appendix  contains  plots  for  one  elevation  and  consists  of 
six  sets  of  three  figures.  Each  set  corresponds  to  a  particular  range: 
34.48,  60.96,  91.44,  152.4,  304.8,  or  609.6  meters  from  ground  zero. 

For  each  range  fluence,  flux,  and  average  flux  are  presented  as  a  func- 
tion  of  time.  Fluence  is  in  units  of  cal -cm  ,  flux  and  average  flux 
are  in  units  of  cal -cm-2-sec  . 
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Figure  85.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  30.48  meters  ground 
range. 
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Figure  86.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  60  degrees  at  30.48  meters  ground 
range. 
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Figure  87.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  60  degrees  at  30.48  meters  ground  range 
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Figure  88.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  60.96  meters  ground 
range. 


163 


1.0&-04 


Figure  89.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  60  degrees  at  60.96  meters  ground 
range. 
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Figure  90.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  60  degrees  at  60.96  meters  ground  range. 
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Figure  91.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  91.44  meters  ground 
range. 


166 


FLUX  (CAL/ CH  /$EC) 


TIME  SEC 


Figure  92.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  60  degrees  at  91.44  meters  ground 
range. 
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Figure  93. 


Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  91.44  meters  ground  range. 
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Figure  94.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  152.4  meters  ground 
range. 
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Figure  95.  Flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  60  degrees  at  152.4  meters  ground 
range. 
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Figure  96.  Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  152.4  meters  ground  range 
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Figure  97.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  304.8  meters  ground 
range. 
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Figure  98.  Flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  304.8  meters 
ground  range. 
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Figure  99.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  60  degrees  at  304.8  meters  ground  range. 
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Figure  100.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  609.6  meters  ground 
range. 
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Figure  101.  Flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  609.6  meters 
ground  range. 
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Figure  102. 


Average  flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  60  degrees  at  609.6  meters  ground  range. 
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APPENDIX  D 


ESTIMATE  OF  THERMAL  RADIATION 
FROM  A  1  KILOTON  NUCLEAR  SURFACE  BURST 
DETECTOR  ELEVATION  90  DEGREES 


This  appendix  contains  plots  for  one  elevation  and  consists  of 

six  sets  of  three  figures.  Each  set  corresponds  to  a  particular  range: 

30.48,  60.96,  91.44,  152.4,  304.8,  or  609.6  meters  from  ground  zero. 

For  each  range  fluence,  flux,  and  average  flux  are  presented  as  a  func- 

-  2 

tion  of  time.  Fluence  is  in  units  of  cal -cm  ,  flux  and  average  flux 
are  in  units  of  cal-cnf 2-sec-1 . 
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Figure  103.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  elevation  90  degrees  at  30.48  meters  ground 
range. 
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Figure  104.  Flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  30,48  meters 
ground  range. 
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Figure  105.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  90  degrees  at  30.48  meters  ground  range. 
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Figure  106.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  60.96  meters  ground 
range. 
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Figure  107.  Flux  from  1  *  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  60.96  meters 
ground  range. 
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Figure  108.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  90  degrees  at  60.96  meters  ground  range. 
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Figure  109.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  91.44  meters  ground 
range. 
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Figure  110.  Flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  91.44  meters 
ground  range. 
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Figure  111.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  90  degrees  at  91.44  meters  ground  range. 
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Figure  112.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  152.4  meters  ground 
range. 
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Figure  113.  Flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  152.4  meters 
ground  range. 
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Figure  114.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  90  degrees  at  152.4  meters  ground  range. 
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Figure  115.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  304.8  meters  ground 
range. 
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Figure  116.  Flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  304.8  meters 
ground  range. 


193 


AVERAGE  FLUX  (CAL/CM7SEC) 


— - 1  I  I  I  i  II  ;  i  r  I  I  I  ll  l  l  i  I  I  I  i  i  i  i  i 

1.0E-03  1.0E-02  l.OE-Ol  l.OE+OO 


TIME  SEC 


Figure  117.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  90  degrees  at  304.8  meters  ground  range. 
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Figure  118.  Fluence  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  609.6  meters  ground 
range. 
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Figure  119.  Flux  from  1  kt  surface  burst  crossing  detector 
with  orientation  90  degrees  at  609.6  meters 
ground  range. 
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Figure  120.  Average  flux  from  1  kt  surface  burst  crossing  detector 

with  orientation  90  degrees  at  609.6  meters  ground  range. 
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BMDSC,  E. 

Williams 

4  cy  ATTN: 

TITL 

ATTN 

BMDSC-HW 

4  cy  ATTN: 

RAAE,  W.  McKechney 

ATTN 

BMDSC-H,  N 

.  Hurst 

ATTN 

BMDSC-NW 

Defense  Technical  Information  Center 

ATTN 

BMDSC-HLA, 

R.  Webb 

12  cy  ATTN: 

DD 

Chief  of  Engineers 
Department  of  the  Army 
ATTN:  DAEN-MCE-D 


Defense  Nuclear  Agency 

Department  of  the  Army 

ATTN 

FCTMD 

ATTN: 

CERL-SOI-L 

ATTN 

FCPR 

ATTN 

FCTMOT 

Deputy  Chief  of  Staff  for 

Rsch  Dev 

ATTN 

FCTMD,  W.  Summa 

Department  of  the  Army 

ATTN 

FCTMOF 

ATTN: 

DAMA-CSM-N 

ATTN 

FCPR,  J.  Hill 

ATTN 

FCTMEI 

Electronics 

Tech  &  Devices 

i  Lab 

ATTN 

FCT 

U.S.  Army  Electronics  R&D 

Command 

ATTN: 

DRSEL 

Field  Command 

ATTN: 

DELET-ER 

Defense  Nuclear  Agency 

ATTN: 

WPNS  Effects 

Section 

ATTN: 

:  FCPRL 
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DEPARTMENT  OF  THE  ARMY  (Continued) 


DEPARTMENT  OF  THE  ARMY  (Continued) 


Harry  Diamond  Laboratories 

Department  of  the  Army 
ATTN:  DELHD-I-TL 
ATTN:  DELHD-N-P,  L.  Belliveau 
ATTN:  DELHD-N-P,  F.  Wimenitz 
ATTN:  DELHD-N-P 
ATTN:  DELHD-N-TS 
ATTN:  DELHD-N-P,  J.  Meszaros 

U.S.  Army  Armament  Material  Readiness  Command 
ATTN:  MA  Library 

U.S.  Army  Armament  Rsch  &  Dev  Command 
ATTN:  DRDAR-LCW 

U.S.  Army  Ballistic  Research  Labs 
ATTN:  DRDAR-BLT,  J.  Keefer 
ATTN:  DRDAR-BLT,  N.  Schuman 
ATTN:  DRDAR-BLT,  W.  Taylor 
ATTN:  DRDAR-TSB-S 

U.S.  Army  Cold  Region  Res  Engr  Lab 
ATTN:  Technical  Director 

U.S.  Army  Conn-Elec  Engrg  Instal  Agency 
ATTN:  Tech  Library 

U.S.  Army  Missile  R&D  Command 
ATTN:  DRDMI-TKP 

U.S.  Army  Mobility  Equip  R&D  Cmd 
ATTN:  DRDME-WC 

U.S.  Army  Nuclear  &  Chemical  Agency 
ATTN:  Library 
ATTN:  Library  for  ATCA-NAW 
ATTN:  Library  for  MONA-SAL 
ATTN :  MONA-WE 

U.S.  Army  Tank  Automotive  R&D  Command 
ATTN:  DRDTA-UL 

U.S.  Army  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TDC ,  R.  Benson 

U.S.  Army  Training  &  Doctrine  Cmd 
ATTN:  ATCD-T 
ATTN:  ATORI-OP 

White  Sands  Missile  Range 

Department  of  the  Army 
ATTN:  STEWS-FE-R 
ATTN:  Library 

U.S.  Army  Communications  Command 

ATTN:  Technical  Reference  Division 
ATTN:  CC-OPS-WR,  R.  Nelson 
ATTN:  CC-OPS-PD 

U.S.  Army  Communications  R&D  Command 
ATTN:  DRDCO-COM  RM,  L.  Dorkin 

U.S.  Army  Edgewood  Arsenal 

ATTN:  SMUEA-TS-TM,  R.  Navin 

U.S.  Army  Engineer  Center 
ATTN :  DT-LRC 


U.S.  Army  Engineer  Div  Huntsville 
ATTN:  HNDED-SR 

U.S.  Army  Engineer  Div  Ohio  River 
ATTN:  ORDAS-L 

U.S.  Army  Engineer  School 
ATTN:  ATZA-DTE-ADM 

U.S.  Army  Engr  Waterways  Exper  Station 
ATTN:  WESSS,  J.  Ballard 
ATTN:  WESSE,  L.  Ingram 
ATTN:  Library 
ATTN:  WESSA,  W.  Flathau 

U.S.  Army  Foreign  Science  &  Tech  Ctr 
ATTN:  DRXST-SD 

U.S.  Army  Materiel  Dev  &  Readiness  Cmd 
ATTN:  DRXAM-TL 
ATTN:  DRCDE-D 

U.S.  Army  Missile  Command 

ATTN:  DRSMI-MSM,  D.  Loney 
ATTN:  Library 

DEPARTMENT  OF  THE  NAVY 

David  Taylor  Naval  Ship  R&D  Ctr 
ATTN:  Code  L42-3 

Naval  Air  Systems  Command 

ATTN:  AIR-360G  JP-2,  J.  Schultz 

Naval  Coastal  Systems  Laboratory 
ATTN:  D.  Sheppard 

Naval  Construction  Battalion  Center 
ATTN:  Code  L51,  $.  Takahashi 
ATTN:  Code  L08A 
ATTN:  Code  L51,  R.  Odell o 

Naval  Electronic  Systems  Command 
ATTN:  PME  117-21 

Naval  Electronics  Systems  Command 
ATTN:  Commander 

Naval  Facilities  Engineering  Command 
ATTN:  Code  04B 

Naval  Material  Command 
ATTN:  MAT  08T-22 

Naval  Research  Laboratory 

ATTN:  Code  5584,  E.  Friebele 
ATTN:  Code  7780 

ATTN:  Code  4709,  W.  Ali 
ATTN:  Code  7550,  J.  Davis 

ATTN:  Code  5584,  G.  Sigel 

ATTN:  Code  6770 

ATTN:  Code  8404,  H.  Pusey 
ATTN:  Code  2627 

Naval  Surface  Weapons  Center 
ATTN:  Code  E21 

ATTN:  Code  X211 
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Naval  Surface  Weapons  Center 
ATTN:  W.  Wi shard 

ATTN:  Tech  Library  &  Info  Services  Branch 

Naval  Weapons  Center 
ATTN:  Code  233 
ATTN:  Code  3263,  J.  Bowen 
ATTN:  Code  266,  C.  Austin 

Naval  Weapons  Evaluation  Facility 
ATTN:  Code  10 
ATTN:  R.  Hughes 

Office  of  Naval  Research 

ATTN:  Code  474,  N.  Perrone 
ATTN:  Code  715 

Office  of  the  Chief  of  Naval  Operations 
ATTN:  OP  654E14 
ATTN:  OP  654C3,  R.  Piacesi 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
ATTN:  NSP-43 

DEPARTMENT  OF  THE  AIR  FORCE 

Aeronautical  Systems  Division 
Air  Force  Systems  Command 
ATTN:  ASD/ENFTV 
ATTN:  AS D/YH-YEF,  Capt  Guice 

Air  Force  Armament  Laboratory 
ATTN:  DLYV,  J.  Collins 

Air  Force  Institute  of  Technology 
ATTN:  Library 

Air  Force  Systems  Command 

ATTN:  Technical  Library 

Air  Force  Weapons  Laboratory 
Air  Force  Systems  Command 
ATTN:  DYV 

ATTN:  NTE,  M.  Plamondon 
ATTN:  NTED,  R.  Matalucci 
ATTN:  NTES-G 
ATTN:  NTES-C ,  R.  Henny 
ATTN:  SUL 
ATTN:  DEX 

Assistant  Chief  of  Staff 
Intelligence 

Department  of  the  Air  Force 
ATTN:  I NT 
ATTN:  IN 

Ballistic  Missile  Office 
Air  Force  Systems  Command 
ATTN:  MNNH 

Deputy  Chief  of  Staff 
Research,  Development,  &  Acq 
Department  of  the  Air  Force 
ATTN:  AFRD 

Foreign  Technology  Division 
Air  Force  Systems  Command 
ATTN:  NIIS  Library 
ATTN:  SDBF,  S.  Spring 


DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN:  RST 
ATTN:  RSS 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN:  SZJ ,  L.  Doan 

Strategic  Air  Command 
Department  of  the  Air  Force 
ATTN:  XPQM 

ATTN:  NRI-STINFO  Library 

DEPARTMENT  OF  ENERGY 

Department  of  Energy 
Nevada  Operations  Office 

ATTN:  Mail  &  Records  for  Tech  Lib 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

Brookhaven  National  Laboratory 
ATTN:  P.  Levy 

Lawrence  Livermore  National  Laboratory 
ATTN:  Tech  Information  Dept 

Los  Alamos  National  Scientific  Laboratory 
ATTN:  M.  Pongratz 
ATTN:  Librarian 
ATTN:  MS  670,  J.  Hopkins 

Oak  Ridge  National  Laboratory 

ATTN:  Central  Research  Library 
ATTN:  Civ  Def  Res  Proj,  Kearny 

Sandia  National  Laboratories 

ATTN:  Library  &  Security  Class  Div 

Sandia  National  Laboratories 
ATTN:  L.  Vortman 
ATTN:  3141 
ATTN:  L.  Anderson 

OTHER  GOVERNMENT 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 

Department  of  Commerce 
National  Bureau  of  Standards 
ATTN:  R.  Levine 

Federal  Emergency  Management  Agency 

ATTN:  Asst  Dir  for  Rsch,  J.  Buchanon 
ATTN:  Hazard  Eva!  &  Vul  Red  Div 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  L.  Selzer 

ATTN:  Technical  Information  Services 

Applied  Theory,  Inc 
2  cy  ATTN:  J.  Trulio 

Artec  Associates,  Inc 
ATTN:  D.  Baum 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


AVCO  Research  &  Systems  Group 
ATTN:  W.  Reinecke 
ATTN:  A.  Pall  one 
ATTN:  Library  A830 

BDM  Corp 

ATTN:  Corporate  Library 
Boeing  Co 

ATTN:  Aerospace  Library 
ATTN:  M/S  85/20,  E.  York 
ATTN:  R.  Holmes 

Boeing  Co 

ATTN:  M/S  42/37,  K.  Friddell 

Charles  Stark  Draper  Lab,  Inc 
ATTN:  Tech  Library 

University  of  Denver 
Colorado  Seminary 
Denver  Research  Institute 

ATTN:  Sec  Officer  for  J.  Wisotski 

Effects  Technology,  Inc 
ATTN:  R.  Wengler 

EG&G  Washington  Analytical  Svcs  Center,  Inc 
ATTN:  Library 

EG&G,  Inc 

ATTN:  P.  Zavaharo 


Electromechanical  Sys  of  New  Mexico,  Inc 
ATTN:  R.  Shunk 

Eric  H.  Wang 

Civil  Engineering  Rsch  Fac 
University  of  New  Mexico 
ATTN:  E.  Wang 
ATTN:  D.  Calhoun 
ATTN:  Technical  Library 
ATTN:  G.  Lane 

Gard,  Inc 

ATTN:  G.  Neidhardt 

General  Electric  Company— TEMPO 
ATTN:  W.  Chan 
ATTN:  DAS  I AC 
ATTN:  J.  Shoutens 

Geocenters,  Inc 

ATTN:  E.  Marram 


Kaman  Sciences  Corp 
ATTN:  D.  Sachs 
ATTN:  Library 

Karagozian  and  Case 

ATTN:  J.  Karagozian 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  TIC-Library 

Los  Alamos  Technical  Associates,  Inc 
ATTN:  P.  Hughes 
ATTN:  J.  Kimmerly 

Lovelace  Biomedical  Envir  Rsch  Inst,  Inc 
ATTN:  Technical  Library 

Lulejian  &  Associates,  Inc 
ATTN:  R.  Daniels 

Management  Science  Associates 
ATTN:  K.  Kaplan 

Merritt  CASES,  Inc 
ATTN:  Library 

Mission  Research  Corp 

ATTN:  Tech  Library 

Mission  Research  Corp-San  Diego 
ATTN:  V.  Van  Lint 

Nathan  M.  Newmark  Consul  Eng  Svcs 
ATTN:  W.  Hall 

National  Academy  of  Sciences 

National  Materials  Advisory  Board 
ATTN:  D.  Groves 

Nichols  Research  Corp,  Inc 
ATTN :  N .  Byrn 

Pacific-Sierra  Research  Corp 
ATTN:  H.  Brode 

Pacifica  Technology 

ATTN:  Tech  Library 

Physics  Applications,  Inc 
ATTN:  F.  Ford 

Physics  International  Co 
ATTN:  J.  Shea 
ATTN:  F.  Sauer 


H-Tech  Labs,  Inc 

ATTN:  B.  Hartenbaum 

I  IT  Research  Institute 

ATTN:  A.  Longinow 
ATTN:  Documents  Library 

Information  Science,  Inc 
ATTN:  W.  Dudziak 

Kaman  AviDyne 

ATTN:  N.  Hobbs 
ATTN:  Library 


R  &  D  Associates 

ATTN:  A.  Kuhl 
ATTN:  J.  Carpenter 
ATTN:  F.  Field 
ATTN:  P.  Rausch 
ATTN:  R.  Port 

ATTN:  Technical  Information  Center 

R  &  D  Associates 

ATTN:  B.  Yoon 

Rockwell  International  Corp 
ATTN:  Library 
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Science  Applications,  Inc 
ATTN:  D.  Hove 
ATTN:  Tech  Library 

Science  Applications,  Inc 
ATTN:  J.  McRary 
ATTN:  W.  Mendes 


Science  Appl ications,  Inc 


ATTN 

R.  Miller 

ATTN 

Tech  Library 

ATTN 

J.  Dishon 

Science  Appl ications,  Inc 

ATTN 

W.  Chadsey 

ATTN 

W.  Layson 

ATTN 

M.  Knasel 

ATTN 

J.  Cockayne 

ATTN 

G.  Binninger 

ATTN 

R.  Sievers 

ATTN 

W.  Moechner 

3  cy  ATTN 

B.  Chambers  III 

Science  Applications,  Inc 

ATTN 

K.  Sites 

Southwest  Research  Institute 

ATTN 

W.  Baker 

SRI  International 

ATTN 

D.  McDaniels 

ATTN 

D.  Keough 

ATTN 

A.  Burns 

ATTN 

G,  Abrahamson 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Systems,  Science  &  Software,  Inc 
ATTN:  Library 
ATTN:  R.  Duff 

TeleDyne  Brown  Engineering 

ATTN:  Technical  Library 
ATTN:  J.  Raven sc raft 


Tetra  Tech, 

Inc 

ATTN: 

Library 

TRW 

Defense 

&  Space  Sys  Group 

ATTN: 

Technical  Information  Ctr 

ATTN: 

N.  Lipner 

ATTN: 

J.  Tambe 

ATTN: 

B.  Sussholtz 

ATTN: 

R.  Eastman 

TRW 

Defense 

&  Space  Sys  Group 

ATTN: 

G.  Hulcher 

ATTN: 

P.  Dai 

foreign  governments 


Admiralty  Surface  Weapons  Establishment 
Ministry  of  Defense 

ATTN:  W.  Del  any 

Atomic  Weapons  Research  Establishment 
2  cy  ATTN:  D.  Sanuels 

British  Defense  Staff 
British  Embassy 
ATTN:  ACOW 
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